Topological liquids and valence cluster states in two-dimensional SU(A^) magnets 
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We study the zero temperature phase diagram of a class of two-dimensional SU(A^) antiferromag- 
nets. These models are characterized by having the same type of SU(A'^) spin placed at each site of 
the lattice, and share the property that, in general, more than two spins must be combined to form 
a singlet. An important motivation to study these systems is that they may be realized naturally 
in Mott insulators of alkaline earth atoms placed on optical lattices; indeed, such Mott insulators 
have already been obtained experimentally, although the temperatures are still high compared to 
the magnetic exchange energy. We study these antiferromagnets in a large- A'^ limit, finding a variety 
of ground states. Some of the models studied here have a valence bond solid ground state, as was 
found in prior studies, yet we find that many others have a richer variety of ground states. Focusing 
on the two-dimensional square lattice, in addition to valence cluster states (which are analogous to 
valence bond states), we find both Abelian and non-Abelian chiral spin liquid ground states, which 
are magnetic counterparts of the fractional quantum Hall effect. We also find a "doubled" chiral spin 
liquid ground state that preserves time reversal symmetry. These results are based on a combination 
of rigorous lower bounds on the large- A^ ground state energy, and a systematic numerical ground 
state search. We conclude by discussing whether experimentally relevant SU(A) antiferromagnets 
- away from the large-A limit - may be chiral spin liquids. 



I. INTRODUCTION 

Since their discovery nearly thirty years ago, fractional 
quantum Hall (FQH) liquids continue to be a rich source 
of novel and exciting physics Ja2 FQH liquids belong to an 
intriguing class of quantum states of matter: they do not 
fall into the conventional classification in terms of bro- 
ken symmetry, electron band structure, and Fermi liquid 
theory, but instead are characterized by the notion of 
topological orderj^i^ In FQH liquids, topological order 
is directly responsible for the celebrated properties of 
fractional charge, fractional and non-Abelian statistics, 
and gapless chiral edge states. While the phenomenon 
of topological order is not limited to FQH systems in 
principle, they remain its only known experimental re- 
alization, recent progress with rotating cold atomic con- 
densates notwithstanding.^ It is thus important to ask in 
which other systems we might find topologically ordered 
states of matter. 

Over the past several years there has been consider- 
able progress identifying model quantum spin systems ex- 
hibiting topological order 1^"— The ground states of these 
models typically have no spontaneously broken symme- 
tries, and are thus referred to as quantum spin liquids; 
these are concrete realizations of Anderson's idea of res- 
onating valence bonds. The models that can be shown 
to exhibit topological order are generally not very re- 
alistic, but many are built from realistic degrees of free- 
dom without any special symmetries, making it clear that 
there is no in-principle obstacle for topological order to 
exist in real quantum spin systems. Despite this progress, 
we know of no candidate materials for a topologically or- 
dered spin liquid. While a few solid state quantum mag- 
nets are quantum spin liquid candidates,—"— all these 
systems seem to have gapless excitations and are thus 
not natural candidates for topological order. 



Here, we discuss a class of spin systems with topo- 
logically ordered spin liquid ground states. While the 
systems we study are not realistic for solid state materi- 
als, they can be realized naturally - without complicated 
engineering of a special Hamiltonian - using fermionic ul- 
tracold alkaline earth atoms (AEA) in optical latticesi^ 
While most ultracold atom experiments to date involve 
alkali atoms, AEA are promising systems to study many- 
body physics, and experiments in this direction are pro- 
gressing rapidlyj^"— An important feature of these sys- 
tems is the presence - without fine-tuning - of a large 
SU(7V) spin- rotation symmetry, where N — 2/+ 1, and / 
is the nuclear spini^i^^ The nuclear spin can be as large 
as / = 9/2 for ®''Sr, so N can be as large as 10. The focus 
of this paper is primarily on the models themselves, and 
not their cold atom realizations; nonetheless, for com- 
pleteness we review in Appendix the realization of the 
spin systems of interest using AEA. Further information 
along these lines can be found in Ref. 2l!. 

The models we study are two-dimensional SU(iV) anti- 
ferromagnets where the SU(A^) representation (i.e. type 
of spin) is the same on every lattice site - the simplest 
case is the iV-dimensional SU(A^) fundamental represen- 
tation. More generally, we consider spins in the SU(iV) 
irreducible representation labeled by a m x ric Young 
tableau with m < N rows and Uc columns (Fig. [T]). We 
will refer to such a representation as the m x ric repre- 
sentation. These models differ crucially from solid state 
SU(2) magnetism, in that, in general, more than two 
spins are required to form a SU(iV) singlet. This means 
that singlets are not two-site valence bonds, but rather 
are multi-site "valence clusters." The cases Uc = I and 
ric = 2 can both be realized as AEA Mott insulators 
(Appendix While a variety of values of m can be 
realized, m = 1 is of greatest interest because it best 
avoids issues of three-body and other losses. For SU(2) 
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FIG. 1. The Young tableau corresponding to the m x Uc ir- 
reducible representation of SU(A''). We consider spin models 
where the spin at each lattice site transforms in this repre- 
sentation. 



spins, TO = 1 and Uc = 25", so that Uc — 1,2 correspond 
to S* = 1/2, 1, respectively. It should be noted that these 
models are distinct from a much-studied class of SU(Af) 
spin models, where (inequivalent) conjugate representa- 
tions occupy the two sublattices of a bipartite latticei^i^i 

In Ref. [3^ together with A. M. Rey, we considered the 
semiclassical limit to = 1 and ri^ — > 00. This is analo- 
gous to the large-S limit for SU(2) spins, and, indeed, 
reduces to it when N = 2. The limit — >■ 00 is bi- 
ased toward magnetically ordered ground states, because 
the spins become classical A'^-component complex vec- 
tors. However, it turns out that the ground state man- 
ifold is in general extensively degenerate (precisely, its 
dimension is proportional to the number of sites in the 
system) - on the square lattice with nearest-neighbor ex- 
change this occurs for > 3, with the degree of exten- 
sive degeneracy growing with N. This situation some- 
times occurs in geometrically frustrated magnets, where 
a common consequence is that magnetic order is strongly 
suppressed, and sometimes even destroyed, by large very 
low-energy fluctuations.— Given that this occurs even 
in a limit which is deliberately biased in favor of mag- 
netic order, non-magnetic ground states are likely for the 
ric — 1,2 cases of greatest interest. (We note that recent 
work has given strong evidence that a magnetically or- 
dered ground state does occur for m — Uc — 1 and A^ = 3 
on the square lattice 1^ For A^ = 4, while prior exact 
diagonalizatioii2i, and variational wavefunction^S studies 
favored a non-magnetic ground state, a very recent study 
employing both projected entangled pair states and ex- 
act diagonalization found evidence in favor of magnetic 
order4S These results are consistent with the expecta- 
tion that non-magnetic ground states are more likely for 
larger values of A^, where the extensive degeneracy in the 
semiclassical limit is larger.) 

This paper is concerned with the ground states of these 
SU(A^) antiferromagnets, in a solvable large- A^ limit suit- 
able for addressing the competition among non-magnetic 
states4ii^ With Rey in Ref. [H, we studied the case 
= 1 on the square lattice in the large- A^ limit, and 
announced a number of results. Here, we study the case 
of arbitrary Uc on general lattices, with a focus on the 
square lattice for ric = 1, 2. We also provide more detail 
on the results already reported in Ref. [ssl 

In the large- A^ limit, A^ and m are taken to infinity. 



while N/m = k and Uc are held fixed. The parameter 
fc, which we choose to be an integer greater than unity, 
plays a very important role in our analysis: k is the min- 
imum number of spins needed to form a SU(A^) singlet. 
Given this physical interpretation of fc, the large- A^ limit 
can thought of as a solvable generalization of the model 
with SU(fc) symmetry, to = 1, and the same fixed Uc- 
Readers primarily interested in the implications of our 
results for real AEA Mott insulators can interpret our 
large-A^ results as a prediction for the ground state of 
these physically realizable m — 1 models. This bold pre- 
diction will need to be tested further in future work; see 
Sec. IVIII for further discussion along these lines. 

In general, the SU(A') singlets are fc-site valence clus- 
ters. Based on the observation that, when fc = 2, 
singlets are 2-site valence bonds, the case k = 2 has 
been studied as a solvable large-A^ generalization of 
SU(2) antiferromagnetism<^^i^"— For the same reason, 
the k > 2 case does not provide a good generalization 
of SU(2) antiferromagnetism. Under very general con- 
ditions in the k — 2 large-A^ limit, the ground state is 
a valence-bond solid (VBS) that spontaneously breaks 
lattice symmetries!^ One of the striking results of this 
paper (and Ref. HI) is that the large- A" ground states are 
much richer in the less-studied case k > 2. 

While SU(A') spin models with the same representa- 
tion on every lattice site have not received extensive at- 
tention (except in the case of self-conjugate representa- 
tions, i.e. fc = 2), there have been several earlier studies. 
While our focus is primarily on two dimensions, we note 
that the one-dimensional chain with to = ric = 1 was 
solved exactly for all A^;^ and the effective field theories 
of it and other chains were also studied. In two di- 
mensions, most work focused on the m, — Uc — 1 model 
with either A^ = 3 or A^ = 4. The former case arises 
as a special point of a 5* = 1 spin model with bilinear 
and biquadratic exchange terms while the latter is 
a highly symmetric point of a S = 1/2 Mott insulator 
with an additional two-fold orbital degeneracyj^Sii^i^i^ 
or a special point of a model with Sp(4) symmetry4^ 
We also note a further very recent study of the A^ = 4, 
TO = = 1 model on the square lattice4S Models on the 
cubic lattice have been studied in high-temperature series 
expansion,^" and a class of exactly solvable models with 
TT-c > 1 was studied in Ref. [Sll Finally, effective models 
of valence cluster degrees of freedom - analogous to more 
familiar quantum dimer models - have been studiedi^>^ 

Returning for a moment to the ultracold atom realiza- 
tion of our models, it should be mentioned that high-spin 
quantum magnets can also be realized using alkali atoms, 
and in that context also have spin symmetry enhanced 
above SU(2)42. However for an A^-component system, the 
symmetry is generically less than SU(A^). For example, 
in a spin-3/2 alkali system, the spin symmetry is ex- 
pected generically to be Sp(4) and not SU(4).— While 
these systems share with SU(2) magnets the property 
that two spins can be combined to form a singlet, they 
are also likely to be fertile ground for the realization of a 



3 



variety of interesting ground states j — In close rela- 
tion to quantum magnetism, half-filled repulsive SU(A^) 
Hubbard models have also been studied in the context of 
ultracold atomsi^i^ Very recently, the repulsive SU(3) 
Hubbard model was studied for arbitrary filling. 

We now summarize our results for the square lattice 
with ric = 1,2. (A graphical summary for ric = 1 can 
be found in the phase diagram of Fig. ^ discussed in 
Sec. lVHI ') Depending on fc, we find valence cluster states 
(VCS) that break lattice symmetries and are formed by 
tiling the lattice with multi-site singlet clusters, and three 
distinct types of topologically ordered spin liquids. Two 
of the spin liquids are chiral spin liquid (CSL) states<^"— 
The CSL is a spin system analog of an FQH state; it 
spontaneously breaks parity and time-reversal symmetry 
(symmetries that are broken explicitly by the magnetic 
field in FQH systems),— supports excitations with frac- 
tional quantum numbers and statistics, and has gapless 
chiral edge states that carry spin. 

We find both an Abelian chiral spin liquid (ACSL) (for 
He = 1 and fc > 5), and a no n- Abelian Chiral spin liquid 
(nACSL) (for = 2 and fc > 6). The ACSL is described 
at low energies by a U(l) n Chern-Simons theory and has 
fractional statistics. The nACSL, on the other hand, is 
described by a U(l)2Ar x SU(2)Ar Chern-Simons theory, 
and supports non-Abelian statistics. The third state we 
find is distinct from these CSL states in that it preserves 
time reversal symmetry. At the mean-field level it ap- 
pears as two copies of the ACSL, with opposite chiralities, 
and we thus dub it a doubled chiral spin liquid (dCSL); 
its low-energy theory is a U(l) mutual Chern-Simons the- 
ory. A more concrete way of describing all these states is 
in terms of Gutzwiller projected trial wavefunctions, as 
described in Sec. IHII The dCSL is found for ric = 2 and 
has the same energy as the nACSL in the N ^ oo limit. 
Presumably 1/iV corrections select one of these states as 
a ground state; we have not computed these, since in 
our view the large- A'^ limit is primarily useful as a tool 
to determine likely ground states of physically realizable 
models, and ultimately the issue of whether the dCSL or 
nACSL (or some other state) is lower in energy will need 
to be determined by directly studying those models (oc- 
curring at finite N). We find VCS states for ric = 1,2 and 
2 < fc < 4, as well as a more complicated inhomogeneous 
ground state when Uc = 2 and fc = 5. These results are 
obtained via a combination of exact lower bounds on the 
large- A^ ground state energy (generalizing the results of 
Ref. |43[ ) , and a systematic numerical search. 

A current interest in states which support excitations 
with non-Abelian statistics is fueled by the expectation 
that they could be used to build a topologically pro- 
tected quantum computer.^*' The simplest non-Abelian 
statistics is described by an SU(2)2 Chern-Simons the- 
ory. It is believed to be realized in the quantum Hall 
effect at the filling fraction 5/2,. ^"^'^^^ as well as in a va- 
riety of setups involving Majorana fermionsi^^— How- 
ever, it is not rich enough to support universal quan- 
tum computations.^^ Non-Abelian statistics described by 



SU(2)Ar Chern-Simons theory for iV > 2 is significantly 
richer and in fact gets richer as N increases. In particu- 
lar, A^ = 3 or A^ > 5 is known to be sufficient for univer- 
sal quantum computations.''^ Some fractional quantum 
Hall states in the first excited Landau level are believed 
to realize these types of non-Abelian statistics for mod- 
erate A^, at least for A^ = We observe that the 
non-Abelian statistics proposed here can be as high as 
SU(2)io (in case of *^Sr), thus it is inherently very rich. 

We note that solvable spin models with CSL ground 
states have been found previously)^""— One of these^^ 
is a generalization of the Kitaev model on a decorated 
honeycomb lattice. The models of Refs. Issl - iSGl i nvolve 
long-range 6-spin interactions. Refs. HO and ISlI found 
that a CSL was the large- A^ ground state of a SU(Af) 
spin model of a different type from those considered here, 
where in addition a four-spin ring exchange term which 
explicitly broke time-reversal invariance was added to the 
Hamiltonain. Very recently, also motivated by ultracold 
alkaline earth atoms, Szirmai et. al. studied the same 
type of spin model discussed in this paper, for ric = 1 
and fc = 6 on the honeycomb lattice, and found a CSL 
ground state in the large- A^ limit 

We now give an outline of our paper. In Sec. |TT1 we de- 
fine a broad class of SU(A^) Heisenberg models in terms 
of slave fermions; this is convenient for understanding the 
large- A^ limit, which is also described in this section. In 
Sec. mil we discuss the properties of the Abelian chiral 
spin liquid, non-Abelian chiral spin liquid, and doubled 
chiral spin liquid, including their wavefunctions and edge 
states. We spend the rest of the paper arguing that these 
states indeed appear in the large A^ limit of the appro- 
priate models. In particular, in Sec. IIVI we discuss the 
solution to the large A^ limit of our models on general 
lattices. We give examples of lattices where the large A^ 
solution can be proven to be a VCS, by generalizing re- 
sults of Ref. 43 to fc > 2. The principal tool of analysis is 
a rigorous lower bound on the large- A^ ground state en- 
ergy, which is saturated by certain VCS states. We also 
give general arguments that VCS are not the only states 
which are possible on generic lattices, and other states, 
including spin liquid states, should naturally appear in 
appropriate cases. In Sec. |V] we specialize to bipartite 
lattices, showing that a stricter lower bound on the en- 
ergy can be obtained in this case (when fc > 2). Finally, 
in Sec. lVII we further specialize to the square lattice. Us- 
ing the rigorous lower bounds, we show that the large- A^ 
ground state is a VCS for fc = 2, 3, 4, for both ric = 1 and 
Tie = 2. Next, employing a numerical analysis we show 
that the large- A^ ground state at ric = 1 on the square 
lattice is the ACSL for 5 < fc < 8. Moreover, at ric = 2 we 
show that nACSL and dCSL are the degenerate ground 
states at fc = 6,7. Closely tied to these results is the 
discussion of Appendix [0 where we discuss the possible 
ground states in the limit of large fc. In particular, for 
TT-c = 1, we show that the ACSL wins over VCS states 
as well as a trial uniform gapless state, giving us ammu- 
nition to conjecture that ACSL is the ground state for 
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all fc > 5. The same analysis goes over to ric = 2 and 
leads us to conjecture that the nACSL and dCSL are 
degenerate ground states for all fc > 6. 

We would not have studied these models if it were not 
for the strong potential to realize them in systems of AEA 
on optical lattices. The paper concludes with a discus- 
sion in Sec. IVIIl focusing on the prospects to find chiral 
spin liquids in those spin models that can be realized 
in cold atom experiments. In particular, we discuss the 
phase diagram in the k-m plane (Fig [51). Finally, we men- 
tion some directions for future study; one such direction 
is to understand how fractional or non-Abelian particles 
may be localized and braided in these systems, with an 
eye toward detection of fractional or non-Abelian statis- 
tics. In Appendix If] we further discuss some ideas in this 
direction, describing how fractional holons (which carry 
conserved atom number but not spin), may be localized 
by applying an external potential. 

In Appendix |X1 we review some aspects of the cold 
atom realizations of these systems. Moreover, starting 
from the Hubbard model describing AEA on an optical 
lattice in the large-?/ limit, we derive the appropriate 
Heisenberg models using degenerate perturbation the- 
ory. Some technical details are given in Appendices [B] [Cl 
and[D] 

II. MODELS AND LARGE-iV LIMIT 

Here we introduce the SU(Af) spin models and con- 
struct the solvable large- iV limit, which allows us to 
address the competition among non-magnetic ground 
states. We shall define the models in terms of the 
fermionic spinon operators f^-aa- Here r labels lat- 
tice sites, a — l,...,iV is the SU(7V) spin index, and 
a — I, . . . ,nc will be called a "color" index. The a index 
transforms in the fundamental representation of SU(A^) 
spin rotations; that is, a global SU(A^) rotations acts by 

fraa Uafifrap, (1) 

where U is an arbitrary SU(iV) matrix. (Here, and 
throughout the paper, summation over repeated indices 
is implied. This does not apply to repeated site labels 
r.) Similarly, the a index transforms in the fundamental 
representation of SU(nc) color rotations. It is important 
to distinguish the spinous from the physical fermions of 
an underlying Hubbard model (as in Appendix \^ ; we 
elaborate on this distinction and its importance below. 

Before defining the Hamiltonian, we must specify the 
type of spin at each lattice site. This is accomplished by 
a pair of local constraints, 

flaJvacc = n^Tn (2) 
fraor^abfrba ~ ^- (^) 

Here, A = 1,. . . — 1 labels the traceless, Hermitian 
ric X Tic matrices T"* that generate infinitesimal SU(nc) 
rotations. The proper interpretation of these constraints 



is that, for each lattice site, we restrict to the subspace 
of the fermion Hilbert space spanned by eigenstates of 
the left-hand sides of Eqs. (|2I3|) . with eigenvalues given 
by the right-hand sides. The first constraint specifies 
a fixed number of fermions on each lattice site, and the 
second constraint dictates that each site is a color singlet. 
The second constraint is omitted when ric = 1. These 
constraints project out the "charge" (conserved number) 
and color degrees of freedom of the spinous, which are not 
physical at the microscopic level but are important for 
understanding the low-energy effective theories obtained 
in the large- iV limit. While the constraint Eq. ^ may 
appear mysterious, in the case ric = 2 it arises naturally 
in the large-C/ limit of the Hubbard model describing one 
type of AEA Mott insulator, as described in Appendix [Xl 

Taken together, the constraints imply that the spin at 
each site transforms in the SU(A^) irreducible represen- 
tation with a m X Tie rectangular Young tableau - this 
is shown in Appendix [B] Since all physical operators 
must commute with these local constraints, which to- 
gether form a U(nc) algebra, in this choice of variables, 
there is a local U(nc) redundancy. This is intimately 
related to the fact that, in the large-A^ limit, the low- 
energy effective theory is a U(nc) gauge theory; we shall 
see this below. This type of slave particle representation 
has been employed before J^i^ii^ 

The SU(A^) spin operators are defined to be 

Sal3{r)=^fLafraP' (4) 
a 

and the Hamiltonian is 

J2 Jrr'S^fi{r)Sp^{r'). (5) 

(r,r') 

Here, the sum is over all pairs of sites {r,r'). The cases 
Tic — 1,1 are realizable with alkaline earth atoms, as dis- 
cussed in Appendix|Xl and most of our analysis is focused 
on these cases. We shall always consider m = N/k, where 
fc > 2 is an integer. The parameter fc, as introduced in 
Sec. m is the minimum number of spins required to form 
a SU(A^) singlet. This model becomes exactly solvable 
in the limit where N and m are taken to infinity, while 
fc and Tie are held fixed. For technical convenience, we 
also write Jrr' = Jrr' /N and hold jTrr' fixed; this cor- 
responds merely to multiplication of the Hamiltonian by 
a constant. The case of greatest experimental interest is 
TO = 1, and the large- limit with fixed k and n,, should 
be thought of as a solvable limit of the model with m = 1 
and N = k; the minimum number of spins required to 
form a singlet is the same as in this model. 

We now describe in detail the large- A^ solution, which 
follows the work of Affleck and MarstoEk^>^, and also 
Read and Sachdevi^ Affleck and Marston studied the 
case Tic — 1 and m = N/2, while Read and Sachdev 
generalized their results to arbitrary ric while still fixing 
TO = A^/2. The formal structure of the large- A^ solution 
is the same as in the earlier works, but, as is discussed in 
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the following sections, the nature of the ground states is 
dramatically different. 

We first consider separately the case nc = f for its 
greater simplicity. The starting point is the imaginary- 
time functional integral for the partition function 



Z = J VfVfVxVX exp ( - S{f, /, x, A)) , 



(6) 



where the action is 



S=[ dry^lfradrfra + 'i>^r{frafra - m) (7) 

Jo ^ 

Jo Jrr' 

Jo I ,^ 

(r,r') 

Here, the fermionic variables fra{T) and fra{T) are the 
usual Grassmann variables. \r{T) is a real Lagrange mul- 
tiplier field that implements the constraint flafra — ™- 
The primed sum in the last two terms is over 

only those bonds {r,r') where Jrr' 7^ 0, and Xrr'{T) is 
a complex field defined on the same set of bonds. Upon 
integrating out x one obtains the Hamiltonian Eq. ([5]), 
which is quartic in fermion operators. We focus on the 
zero-temperature limit /3 — ?> cxd. 

We can formally integrate out the fermions and obtain 
an effective action 



Scs{xA) = / 

Jo , „, 



(r ,r') 

7VTrlnQ(x,A), 



'''^—\Xrr'f - im 

Urr' 



rl3 
-'0 



(8) 



where Q is the quadratic form characterizing the 
fermionic part of the action Eq. ([T]). Since m = N/k, SeS 
has a prefactor of N and no other A'^-dependence, imply- 
ing that when N ^ 00 the functional integral over x and 
A can be done exactly using the saddle point approxima- 
tion. We therefore replace x £^nd A by non-fluctuating 
fields 



Xrr' ^ Xrr' 
A^^ y Tj ^r , 



(9) 
(10) 



which are substituted into Eq. (O to obtain a theory of 
non-interacting fermions subject to the mean-field Hamil- 
tonian 



n 



MFT 



{r.r') 



N 

\Jrv' 



(11) 



^ (Xrr'/Ja/r'a +H-C- 



The imaginary saddle point A^ — > i[ir is needed for the 
mean-field Hamiltonian to be Hermitian. We emphasize 
that despite the appearance of the term "mean-field," 



and the appearance of mean-field equations and a mean- 
field Hamiltonian, we are not making any sort of mean- 
field approximation. That is, in the N 00 limit, the 
specific mean-field decoupling we consider - and only this 
decoupling - becomes exact. The results we present are 
thus exact for the Heisenberg spin model in the N 00 
limit. 

In order for this to be a legitimate saddle point, we 
must satisfy the extremum condition 



-S, 



off 



5\r 



5eff 



0. (12) 



In the low-temperature limit, the ground state energy 
Emft of Hmft satisfies S^eix ^ A -J> i^i) = /3£^mft, 
so satisfying Eq. ([T^ is equivalent to extremizing the 
ground state energy. The saddle point equations of 
Eq. are equivalent to the more convenient expres- 
sions 

_ Jr't 

Xrr' - 



N 



^ ifrafra)- 



(13) 
(14) 



Here, the expectation values are calculated using the non- 
interacting Hamiltonian iJMFT- 

Now that we have discussed the simpler case ric = 1, 
we discuss the general case > 1, describing only those 
aspects that differ from nc = 1, and making some defi- 
nitions that will be useful later on. The principal differ- 
ence is that now the fields x and A are Uc x Uc matrices: 
Xrr' (■'') is a general complex Uc x Uc matrix, and A°^(t) is 
a Uc X He Hermitian matrix. The partition function is of 
the same form as Eq. ([5]), where the integration over A is 
understood to be over the restricted space of Hermitian 
matrices. The action is now 

S — dr Ifraadrfraa i{\^^ fraafrba ^ mtr(A^)) 

+ / dr (15) 

+ / dr V'[Xrr'/™a/r';,a+H.C.]. 

Jo , ,s 
[r.r') 

The traces in this expression are in the color space. The 
field At- is again a Lagrange multiplier, now implementing 
both the constraints of Eq. © and ([3]). Again, integrat- 
ing out Xrr' we obtain the Hamiltonian Eq. ([S]). 
The saddle point values of the fields take the form 



/\.rr 



X 



ah 



A-rv 



(16) 
(17) 



where fir is a Hermitian matrix. The mean-field Hamil- 
tonian is 

El iV ,_t 
^ tr(x;, 

, Urr' 
[r.r' ) 



n 



.^tr(Ai.) (18) 
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where 

{r.r') 
r 

Here, we have defined the color density 



= fT f 

J raaJ rba ' 



(19) 
(20) 

(21) 



Note that we can also write T-Ly 
saddle point equations are now 



Arr' 



\Jr\ 



ra S' 



ab 



N 



r'boc^ ra 



J2r tr(^T^nr)- The 

(22) 
(23) 



When analyzing the mean-field Hamiltonian we shall al- 
ways work in the canonical ensemble for the conserved 
fermion number. 

This discussion shows that finding the ground state 
in the large- -/V limit reduces to finding the saddle point 
with lowest energy Emft- In general this task, while a 
great deal simpler than finding the ground state of the 
original quantum problem, is still nontrivial. We shall 
make progress below using a combination of exact lower 
bounds on i^MFT, and numerical methods to search for 
ground states. 

As mentioned above, it is important to recognize that 
the fraa spinon operators are not the same as the phys- 
ical fermions of the Hubbard model discussed in Ap- 
pendix |X1 We are describing a spin model, and there 
are only SU(A'') spin degrees of freedom. In addition to 
spin degrees of freedom the physical alkaline earth atom 
fermions have degrees of freedom associated with their 
conserved number, as well as with their -^Sq and ^Pq elec- 
tronic states. These degrees of freedom are not present 
in the model; this is appropriate for a low-energy de- 
scription of the Mott insulating states we are describing, 
where excitations associated with these degrees of free- 
dom are gapped. (To describe such gapped excitations, 
one must return to the original Hubbard model.) 

The difference between the spinous and the physical 
fermions is manifest when we consider the fluctuations 
about a mean-field saddle point. This allows us to con- 
struct a low-energy effective theory, which goes beyond 
mean- field theory for a given saddle point. The spinous 
in this effective theory should not be interpreted as a mi- 
croscopic representation of the spins, but as low-energy 
effective degrees of freedom. As we shall see below, the 
spinous are minimally coupled to a fluctuating U(nc) 
gauge field. On the other hand, the physical fermions 
of the underlying Hubbard model are uncharged under 
this U(nc) gauge field and do not couple to it directly. 

We consider fluctuations of the form 



\abf\ ^..<^b I ^ab ( ^\ 



Y°-'', 
A,rr 



ir) 



(24) 
(25) 



where a^r' (t) is a ric x ric Hermitian matrix, so that e*"''''' 
is unitary. While other fluctuations are typically trivially 
massive {e.g. amplitude fluctuations in x), these fluctu- 
ations take the form of a U(ric) gauge field minimally 
coupled to the spinous. Specifically, a^r and arr' form 
the time and space components, respectively, of the fiuc- 
tuating U(nc) vector potential. Gauge fiuctuations can 
and do dramatically modify the properties of the mean- 
field state, and therefore should in general not be ne- 
glected. For example, if the gauge field is in a confining 
phase, then the spinous will not be good quasiparticle 
excitations, as a naive mean-field analysis would suggest 
- this indeed occurs in the VCS ground states. On the 
other hand, in the CSL and dCSL phases, Chern-Simons 
terms for the gauge field are present; this not only pre- 
vents spinon confinement, it converts the spinous from 
fermions into anyons. 



III. PROPERTIES OF TOPOLOGICAL LIQUID 
GROUND STATES 



Anticipating the results on energetics discussed below, 
in this section we discuss the properties of the three topo- 
logical liquid ground states on the square lattice. Since 
the main focus of this paper is on energetics, we shall 
content ourselves primarily with deriving low-energy ef- 
fective field theories for each state, and shall not discuss 
the resulting properties in detail. 



A. Abelian chiral spin liquid 

The Abelian chiral spin liquid (ACSL) occurs for Uc 
1, and corresponds to a mean- field saddle point 



Xrr' = xe ' 
^ir = 0, 



(26) 
(27) 



where x is real and positive, and a^^/ is chosen so that 
27r/fc magnetic fiux pierces each plaquette of the square 
lattice. The band structure consists of k bands, of which 
the lowest is full and the others are empty, resulting in a 
Hall conductance (for the mean-field fermions) of axy ~ 
N. 

To understand the properties of this state it is nec- 
essary to go beyond mean-field level, and couple the 
fermions to the fluctuating U(l) gauge field. However 
some properties can already be understood at mean-field 
level. In particular, we see that parity {i.e. reflection) 
and time reversal symmetries are spontaneously broken, 
while the other symmetries of the square lattice (as well 
as SU(A^) spin rotation) are preserved. To see this, it 
is important to recall that in a slave-particle gauge the- 
ory such as this one, symmetry operations act projec- 
tively on the fermions.** For example, if S* : r ^ S{r) 
is a space-group operation, then acting on a fermion it 
may be supplemented by a general space-dependent U(l) 
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gauge transformation: 

S ■ frc 



S{r)a ■ 



(28) 



An operation S" is a symmetry if and only if it is pos- 
sible to find a gauge transformation such that the 
above transformation leaves the mean-field Hamiltonian 
invariant. For the CSL saddle point, this is nothing but 
the familiar magnetic translation group (expanded to in- 
clude all symmetries, not only translations). Because re- 
flections and time reversal both change the sign of the 
gauge-invariant magnetic flux through each plaquette, 
they are spontaneously broken in the ACSL. Other oper- 
ations leave the flux invariant and are indeed symmetries 
of the ACSL. 

To go beyond mean- field theory, we couple the fermions 
to the fluctuating U(l) gauge field. Since the fermions 
are gapped we can integrate them out, resulting in the 
following imaginary-time continuum effective action for 
the gauge field: 



S 



iN 
An 



1 
2? 



(29) 

This is simply Maxwell-Chern-Simons theory. The co- 
efficient of the Chern-Simons term is determined by 
(Jxy — N, while the coefficient of the Maxwell term is 
non-universal (in the large-A^ limit, it is determined by 
details of the fermion band structure). Various proper- 
ties of the ACSL can be derived from this effective action 
- notably, it implies that the fermions are converted via 
fiux attachment into anyons with statistics angle ndnr/N. 

A different - and particularly concrete - route beyond 
mean-field theory is to construct a wavefunction for the 
ACSL. One starts with the ground state of the mean- 
field Hamiltonian, which for the ACSL is simply an inte- 
ger quantum Hall state with the lowest (lattice) Landau 
level filled. One applies the Gutzwiller projection oper- 
ator P, which simply projects onto the subspace with 
exactly m fermions on every lattice site. By construc- 
tion, — T'l'ipo) satisfies the local constraint Eq. ([2]) 
and is thus a legitimate wavefunction for the spin model. 
Such Gutzwiller projected wavefunctions have been stud- 
ied and discussed in a variety of contexts (for a few exam- 
ples, see Refs. Issl - foil ). and properties of such wavefunc- 
tions can be computed numerically using a Monte Carlo 
techniquci^ It is reasonable to expect that 1-0) should 
correctly capture the properties of the corresponding low- 
energy effective gauge theory; for example, it has been 
shown that a class of projected wavefunctions associated 
with an effective Z2 gauge theory capture the expected 
Z2 topological orderi^iS^ However, this expectation will 
need to be tested by future detailed studies of the wave- 
function. In the present case, the projected wavefunction 
may be a useful tool for future microscopic analysis away 
from the large-A^ limit, in particular to help assess the 
prospects for ACSL in physically realizable models. 

Another important property of the ACSL is the pres- 
ence of gapless chiral edge states, which are described by 



a chiral SU(iV)i Wess-Zumino-Witten (WZW) model. A 
simple argument for this can be given following Ref. lo^ : 
Rather than consider the low-energy effective field the- 
ory of fermions coupled to a gauge field, we consider the 
projected wavefunction described above. Before projec- 
tion, the edge mode consists simply of N chiral fermions. 
Using non-Abelian bosonization, the edge theory can be 
cast as two decoupled theories: a chiral SU(A^)i Wess- 
Zumino-Witten (WZW) model, and a chiral U(l) Lut- 
tinger liquid i^^'^^ This is an instance of spin-charge sepa- 
ration, where the SU(7V) spin degrees of freedom are as- 
sociated with the SU(A^)i WZW model, and the fermion 
"charge" with the Luttinger liquid. Upon projection, the 
"charge" degrees of freedom are removed, and hence so 
is the U(l) Luttinger liquid, while the spin degrees of 
freedom and the SU(A^)i WZW model survive. 

Finally, we mention an alternate route to construct a 
low-energy effective theory for the ACSL that does not 
require integrating out the fermions. This approach is 
based on the Chern-Simons effective theory for Abelian 
quantum Hall statesi^ In this approach, one pays the 
price that the SU(iV) symmetry is broken down to 
U(l)^^^, but this is not expected to affect any topo- 
logical properties of the state. Before coupling to the 
gauge field, each spin species of fermion is in an integer 
quantum Hall state, and the current of the fermions of 
spin a (where a — 1, ... , N) can be represented in terms 
of a U(l) gauge field: 



(30) 



The corresponding integer quantum Hall state is cap- 
tured by a Chern-Simons term for which gives the 
following contribution to the real-time Lagrangian: 



Ca - —efj_^xb"d^b°. 



(31) 



Moreover, the coupling of the fermions of spin a to the 
gauge field is simply given by 



(32) 



Finally, the gauge field has no bare Chern-Simons 
term - it is a Lagrange-multiplier field whose role is to 
make the total U(l) fermion current vanish. (The Chern- 
Simons term derived above for came from integrating 
out the fermions.) Combining the above results, we have 
the Lagrangian in X-matrix form. 



1 



C = —KijAie^^xd.Ai 



An 



(33) 



where I ^ I, . . . , N + I, A]^ = a^,, A' = 6^-1 for / > 1, 
and the (TV + f ) X (iV + 1) A"-matrix is 



K = 



I' 

1 1 



(34) 
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Here, Inxn is the N x N identity matrix, and = 
(1, . . . , 1) is a iV-element vector. Following Ref.lOTl both 
bulk and edge topological properties can be deduced 
from this effective theory. We note that the if-matrix 
has N positive eigenvalues and one negative eigenvalue, 
and thus gives rise to N co-propagating edge modes and 
one counter-propagating mode. The counter-propagating 
mode, and one of the co-propagating modes, are singlets 
under U(l)^~^ spin rotations, and these singlet modes 
generically are expected to acquire a gap, leaving — 1 
gapless co-propagating modes - this is nothing but the 
free boson description of the SU(A^)i chiral WZW model. 



B. Non-Abelian chiral spin liquid 

The non-Abelian chiral spin liquid (nACSL) occurs for 
Uc — 2, and corresponds to a mean- field saddle point 



ab 



srab 



Xe'"-'-'()"" (35) 
0, (36) 

where x is I'eal and positive, and again a^^, is chosen 
so that 2Tr/k magnetic flux pierces each plaquette of the 
square lattice. This state has a U(2) — U(l) x SU(2) 
gauge structure, and upon going beyond mean-field the- 
ory the fermions are coupled to a U(2) gauge field. The 
background magnetic flux is a U(l) flux - the background 
SU(2) flux is zero. The band structure can be thought of 
as k 2A^-fold degenerate bands, where the lowest band is 
filled and all others are empty. The mean-field fermions 
have a Hall conductance a^y = 27V. As above, parity 
and time reversal are spontaneously broken, while other 
symmetries are preserved. 

In the large- limit, the ground state energy of the 
nACSL is precisely twice that of the ACSL. This occurs 
because, at the mean-field level, the nACSL is simply 
two decoupled copies of the ric = 1 ACSL, each with the 
same magnetic flux. 

Upon integrating out the fermions, we obtain the fol- 
lowing action: 



S = ^-^ / drcfr e^^xaij.duax 



2Ni 
An 
iN 
47r 



^ / dTcPre^i,xti 



2i 



ai^d,,ax ~ —afj,a„a\ .(37) 



Here is the U(l) gauge field, — Q^^cr' is the 

SU(2) gauge field (tJ* are the usual Pauli matrices), and 
we omitted the Maxwell terms that are also present. The 
second term is the level-iV Chern-Simons term for the 
SU(2) gauge field, which gives rise to the non-Abelian 
statistics of the nACSL. 

As above for the ACSL, one can construct a wavefunc- 
tion for the nACSL. One proceeds as above, but now 
must apply a projection operator to enforce both the lo- 
cal constraints Eq. ^ and Eq. (jS]). 

The chiral edge states of the nACSL can be under- 
stood in terms of an argument very similar to that given 



above for the ACSL. In mean-field theory, there are 2N 
chiral fermions on the edge of the system. Following 
Affleck,— this free fermion theory can be bosonized to a 
chiral SU(iV)2 WZW model (carrying spin excitations), 
a chiral SU(2)7v WZW model (carrying color), and a chi- 
ral U(l) Luttinger liquid. Now the projection removes 
both the "charge" and color degrees of freedom of the 
fermions, leaving only the chiral SU(7V)2 WZW model. 



C. Doubled chiral spin liquid 

The doubled chiral spin liquid (dCSL) occurs for ric 
2, and corresponds to a mean- field saddle point 







= 0, 



(38) 
(39) 



where x and a^^, are as above. In contrast to the nACSL, 
there is now a SU(2) background magnetic fiux, but no 
U(l) fiux. Following the reasoning of Ref. the pres- 
ence of the nontrivial SU(2) fiux breaks the SU(2) gauge 
structure down to U(l). More precisely, the and 
components of the SU(2) gauge field acquire a mass due 
to the presence of the flux, while the component is 
unaffected. Therefore, for the purposes of understanding 
the low-energy physics, we can drop the and com- 
ponents of the SU(2) gauge field, and consider a theory of 
fermions coupled to the two U(l) gauge fields and a^. 
It should be noted that the special role of a^, as com- 
pared to aj^ and a^, is determined by the choice of gauge 
made in writing Eq. ((38)) - a global SU(2) gauge trans- 
formation can be made to select any desired preferred 
axis. 

In the large- TV limit, the ground state energy of the 
dCSL is again precisely twice that of the ACSL, because 
again, at the mean-field level, the dCSL is two decoupled 
copies of the Uc = 1 ACSL, but now with opposite mag- 
netic fiuxes. This means that in the N ^ oo limit the 
dCSL and nACSL have exactly the same energy. This 
degeneracy is expected to be lifted by 1/iV corrections 
that can in principle be computed; this is left for future 
work. 

The dCSL actually respects time reversal symmetry, 
which is implemented by the operation 



r-.fr 



{icr'^)abfrba- 



(40) 



(This operation can be supplemented as well with a 
SU(A^) rotation, but due to the SU(Ar) symmetry this 
is not essential.) The crucial point is that the gauge- 
rotation in the color space compensates for the fact that 
complex conjugation reverses the fiux. Reflection sym- 
metry TZ : r ^ r' , where r' — (—rx,ry), is similarly 
preserved, and 



■ fraa ^ {^'^ )abfr'bo 



(41) 
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The other symmetries (lattice translations and rotations, 
and SU(A^) spin rotations) are preserved in the dCSL as 
they are in the above two states. The dCSL therefore 
does not spontaneously break any symmetries, in con- 
trast to the ACSL and nACSL. 

Upon integrating out the fcrmions, we obtain the fol- 
lowing mutual Chern-Simons action: 



(42) 



Here we have again omitted the Maxwell terms that will 
also be present; the mutual Chern-Simons term fully gaps 
out both gauge fields, and the Maxwell terms play only 
the quantitative role of setting the scale of the gap to 
gauge field excitations. It should be noted that similar 
spin liquid states, but with an additional non-Abelian 
gauge structure, were considered in Ref . [ill. (There, how- 
ever, the analog of the gauge field was incorrectly 
dropped, and therefore a U(l) mutual Chern-Simons 
term was missed.) It can be seen that this term also 
converts the mean-field fermionic excitations into anyons 
with statistics angle tt^tt/N, which can occur in a time- 
reversal invariant fashion due to the color index. We note 
that the same procedure described for the nACSL can be 
applied here to produce a wavefunction for the dCSL. 

Because the dCSL respects time reversal symmetry, it 
lacks chiral edge states. However, it is interesting to note 
that - when N is odd - the edge states are protected at 
the mean-field level, because the mean-field Hamiltonian 
has a nontrivial Z2 topological invarianlrlfiQ, for odd N. It 
is therefore conceivable that topologically protected edge 
states could survive coupling of the mean-field fermions 
to the fiuctuating gauge fields, and it would be interest- 
ing to study this question. Presumably such protection, 
if it occurs, would only hold if one assumes that no spon- 
taneous breaking of time-reversal symmetry occurs at the 
edge. 

Finally, we can also construct a if-matrix Lagrangian 
for the dCSL as above for the ACSL. We let Aj^ = and 
= a^. Next, for / = 3, . . . , + 2, represents the 
current of fermions with a = 1 and spin a = I — 2 (as in 
Eq.[30l), while for / = N+3, 2N+2, A^^ represents the 
current of fermions with a = 2 and spin a — I — [N + 2). 
Following essentially the same reasoning as in Sec. IIII Al 
we have the (2A^ -f 2) x (2A^ + 2) AT-matrix 



K ^ 



/O 

1^ 



OtVxW 



(43) 



where 



ATxJV 



is the N X N matrix of zeros. 



IV. LARGE-iV LIMIT: GENERAL LATTICES 

In Ref. I43!, Rokhsar derived an exact lower bound 
on the large- A^ ground state energy E'mft, for the case 



ric = 1 and k = 2. He further showed that this bound 
is saturated by a VBS state under conditions that are 
satisfied for the great majority of lattices one encoun- 
ters. More precisely, let us define J'^ax to be the largest 
of the exchange couplings J7rr'- (Note that we do not 
restrict to only nearest- neighbor exchange.) Following 
Rokhsar we say that a lattice is dimerizable with respect 
to jTrnax whcn it is possible to partition the lattice into 
2-site dimers, such that the two sites (r, r') in each dimer 
have J^rr' — i/max- Each lattice site must belong to pre- 
cisely one dimer. For a fixed partition into dimers, let 
the set of bonds (r, r') that connect the two sites of a 
dimer be B. Rokhsar considered the VBS saddle point 
defined by 



Xrr' 
Xrr' 



= 



(r,r') e B 
(r,r')^S' 



(44) 



and showed that it is a ground state (its energy satu- 
rates the bound on E'mft)- Except in the case of disor- 
dered systems lacking translation symmetry, most famil- 
iar lattices (and associated sets of JTrr') are dimerizable 
with respect to jTmax-^'^ Therefore, when k = 2, one has 
to consider a relatively unusual lattice to find anything 
other than a VBS ground state in the large- A^ limit. (See 
Ref. [4^ for an example of a lattice that is not dimerizable 
with respect to jTmax-) 

Here, we generalize Rokhsar's bound to the case of 
arbitrary k and ric (Sec. IIV Ap . In Sec. IIVBI we derive 
necessary and sufficient conditions to saturate the bound. 
Next, in Sec. IIV CI we show that the analog of Rokhsar's 
VBS saddle point is a fc-simplex VCS state, where the 
lattice is decomposed into fc-site simplices (fc-simplices 
for short), in which every site is connected to the other 
k — 1 sites by an exchange coupling J7rr' — JTmax- As 
soon as fc > 2, many lattices cannot be decomposed into 
fc-simplices, and for fc > 5 we show that no lattice can 
be decomposed into fc-simplices without fine-tuning of 
the exchange couplings. Therefore it becomes more and 
more difficult to saturate the bound as fc increases. 



A. Derivation of the bound 

Our starting point is the mean-field Hamiltonian 
"Hmft for general ric [Eq. ((HI)], where ^-nd fi'^'' are 
chosen to satisfy the saddle-point equations Eqs. (|22l23p . 
The bound is derived in two steps: first we will show 
-E'mft > ^MFT (defined below), then we will show 
-E'mft — -Abound- The first step was actually omitted 
in Ref. '43'. While this step should not be omitted even 
in the social case considered there, none of the results 
of Ref. |4j are affected by this omission. 

Recalling the definitions of "Hmft in Eq. ([TS]) and Hk 
in Eq. ([T9)). we begin by defining 



^MFT " 



N 



(45) 
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where Xrr' same as in 'Hmft- That is, we obtain 

^MFT by starting with Hmft and setting to zero. 
The ground state energy of "Hmft ^mft- Note that, 
in general, the ground state of H'mft ^il^ ^'^^ satisfy the 
saddle point equations. 

Now, Emft — ('Hmft), where the expectation value is 
taken using the ground state of "Hmft- Using Eq. (f^. 
we note that (Hv) = —n^■J2r^^i^^■>•)^ this cancels the 
second term in "Hmft, so we have 



E 



MFT 



(r,r') 



Jr 



TV /\^T"r 



{H 



(46) 



Letting Ek be the ground state energy of %k, we have 
("Ha') > Ek^ and so 



-Emft > ^mft ~ ^ 



E 



K- 



(47) 



The bound originates from the pair of inequahties 



(54) 
(55) 



which just express the fact that variance is positive. 
These inequalities are saturated (become equalities) if 
and only if is constant over each of the sets C and U. 
Multiplying Eq. ^ by m/N, Eq. ([55]) by {N - m)/N, 
and adding the two, we have 



{N - m) 



Using Eq. (fSSj) and the fact that [e]^ < 0, we have 



\e]r > 



N ~m 



(56) 



(57) 



Now, Ek = ncmNs[e]c, so we have shown 



This is the first of the two desired inequalities. 

We shall now deal with Hf^prp and E^^rj,, and estab- 
lish a lower bound on E^^j,. To do this, we generalize 
Rokhsar's argument^'^ to the case of general m and ric- 
Let Ng be the number of sites of our lattice. We label 
the single-particle energy levels of Hk by an index q; the 
energies are eq. Hk is specified by the UcNNg x ricNNs 
Hermitian matrix 



(48) 



where 



Ek > —UcmNs 



N -m 



(58) 



The next step is to get a simple expression for [e^]. We 
have 



UcNN: 

2 



1 \ ^ 2 



1 



tr[i7|] 



Therefore we have the inequahty 



(59) 



Xr'r — Xrr' ■ 



(49) 



Because this is traceless (all diagonal entries are zero), 
we have 



(50) 



The ground state of "Ha (and hence of H^^prp) is obtained 
by filling the lowest UcmNs energy levels with fermions. 
We call the set of such energy levels C. The other nc{N — 
m)Ns levels, which we denote by the set U, are empty. 

It will be useful to define averages over the sets of levels 
C and U : 



UcmNs 



— E' 



(51) 
(52) 



We also denote the average of over the two sets by 
[e'^]c and [e^]iY, and the average of over all states is 
written [e^]. Equation ([50|) implies 



l^h = - 



N -m 



(53) 



i?WT>^ E'E 

(r,r') a.b 



\ ab 12 
lAr-T"' I 



(60) 



N — m 



m 



\ 



(r,r') a,b 



The next step is to minimize this lower bound, which 
we do by taking the derivative of the right-hand side of 
Eq. ([50)1 with respect to Ix^r'l and setting it to zero: 



= 2A^ 



I Xrr*' I 
\Jrr' 



2m^[N ~m)lm\xfr'\ 



\y.cd |2 

/_^lr''r"') L-^c.d \/^r"r"' 



(61) 



For a given bond (r, r'), this equation implies that either 
IXrM^Oforall a,fe, or 



TLcNs 



E Ei^"v"'i^- 



(r" ,r"') c.d 



m{N-m) 2 

^2 Jrr'- 



(62) 



Now, the left-hand side of Eq. (l62t is independent of the 
bond (r,r'), and so we must have 



ricNs 



E EixX^"i^ = 



2 _ MN - to) „2 



(r",r"') c.d 



7V2 



(63) 
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for some constant Moreover, this implies that, for a 
given bond, unless Jrr' = J*^ then we must have x^r' ~ 
(for all a, 6). Therefore 



E'E^-|EEIx:^M' (64) 

(r,r ) a,o [r.r') a,b 



HcNs m{N — m) 



2 7V2 

Putting these results into Eq. ((60|) . we have 

, "c^s m(iV - m) ^ 

^MFT — >-^*- 



J*. (65) 



2 TV 



(66) 



The global minimum is clearly achieved when = jTnmx, 
the largest of the Jrr' ■ Therefore 



, ficNa m(N — m) 

-C'MFT ^ ^MFT — 7i ^T? Jmax- 



2 N 



Putting m = N/k we have 



(k-l) 

EmFT > -UcNNs Jmax, 



(67) 



(68) 



which reduces to Rokhsar's result when k — 2 and rir — 1. 



Necessary and sufficient conditions to saturate 
the bound 



B. 



Here, we show that the bound Eq. (|6S1) is saturated if 
and only if the following two conditions hold: (1) £q is 
constant over each of the sets £ and U. That is, all the 
filled states have the same energy, and all empty states 
have the same energy. (2) The color density calcu- 



lated using Hk satisfies the condition 



tidirnr) — 0. 



This color density is defined by 



{Ur )K, 



(69) 



(70) 



where the expectation value is taken using the ground 
state of T-Ik- Note that hr in general does not satisfy the 
constraint Eq. 1^^. These conditions for saturation are 
very restrictive, as we discuss below. 

There are two separate inequalities that must both be 
turned into equalities for the bound to be saturated. The 
first is -Emft — -Abound, and the second is -Emft > ^mft- 
Saturation of the first and second inequalities leads to 
conditions (1) and (2) above, respectively. It is trivial to 
show that the first inequality is saturated if and only if 
eq is constant over each of the sets £ and U. 

We now show that condition (2) is equivalent to sat- 
uration of the second inequality. It will be useful to de- 
fine a continuous family of Hamiltonians parametrized by 
ae [0,1]: 



"Ha — T~L 



K 



an 



(71) 



This interpolates between "Hk at a = and "Hk + , 
the fermionic part of 'Hmft, at a = 1. The ground state 
oi Ha with energy Ea is denoted by IV'a)- Because we 
work in the canonical ensemble for the fermion number, 
we are free to make a constant shift /i^^ — > /i"'' + c6'^^ so 
that X^T-trl/^r) = (note that this shift does not change 
Emft)- With this choice for /i^, we have 



(V'll^vlV'i) =0. 



(72) 



We also have ^^mft — -S'mft = E!i — Eq. In particular, 
-Emft = ^^mft if ^^nd only if i^o = -^i . 

The variational principle implies ("(Aal^Q!' IV-'q) ^ Ea'. 
The left-hand side of this inequality can be written 

(V'a|H„'|V'a) =Ea + {a' - a){^a\nv\i'a)- (73) 



We have thus shown 



Ea + {a - a){'lJja\'Hv\^l'a) > Ea 



(74) 



If we put a = 1, this gives Ei > Ea- On the other hand, 
putting a' — 1 gives instead 

Ea + {l-a) {iPa |Hy I Va) > ^1 ■ (75) 

Combining these together, 

Ea + {l-a){'^a\Hv\ii^a)>Ei>Ea, (76) 
which immediately implies 

>0. (77) 

A special case of Eq. ([76|) is 

Eo + {Wv\^o)>Ei>Ea. (78) 

From this it follows that if (V'ol^y IV'o) — 0, then Ei = 
Eq. Now suppose the converse, i.e. suppose Eq = Ei. 
Note that first-order perturbation theory gives us 



and so 



dE 

= {'4'a\'Hv\tpa), 



Ei-Eq= / daiipalHvHa)- 



(79) 



(80) 



By assumption this integral is equal to zero. Since 
the integrand is nonnegative, then we must have 
{ipa\'Hv\''Pa) = 0, and in particular for a = 0. 

Therefore we have shown that Ei = Eq if and only if 
(V'ol'Wy iV'o) = 0, and hence -Emft — Ej^p-p if and only if 
(i/'ol'HvlV'o) = 0. Since (V'ol^vlV'o) = I]rtr(^rnr), we 
have established condition (2) as desired. 

Both conditions derived above for saturation of the 
bound are highly restrictive. Condition (1) dictates that 
there be only two energies in the spectrum; we should 
expect this to occur only when is such that the lat- 
tice is broken into clusters, so that the spectrum consists 
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of perfectly flat bands. Condition (2) is also very restric- 
tive. An easy way to satisfy (2) is simply to have a saddle 
point where fir — 0. Suppose instead that fir 7^ 0, and so 
generically we should expect that fir is non-uniform and 
does not satisfy Eq. ((23)) . It is useful to imagine starting 
from a = and turning on by increasing a. The /ir 
need to be chosen to "even out" the color density, so that 
it satisfies Eq. once a = 1. Naively, a choice of fir 
accomplishing this will cost energy at each lattice site; 
that is, 

(?/'o|[-tr(Airri.)]IV'o) > 0. (81) 
This would imply 

^tlifirflr) < 0, (82) 
r 

which is in conflict with condition (2). This discussion 
indicates that satisfying condition (2) when /i^ 7^ is 
unlikely. 

C. Saturation of the bound and fc-simplex VCS 
states 

The necessary and sufficient conditions derived above 
still leave open the questions of what kind of saddle points 
saturate the bound, and whether saturation is possible 
for a given lattice and set of exchange couplings Jrr' ■ 
Saturation is not always possible - for example, on any 
bipartite lattice with fc > 2, the stricter bound derived 
in Sec. |V] shows that saturation of Eq. (j68| is impossi- 
ble. Here, we will show that, when they exist, fc-simplex 
VCS states saturate the bound and are thus the analogs 
of the VBS states for k = 2. In striking contrast to VBS 
states, many commonly encountered lattices do not ad- 
mit any /c-simplex VCS states for k > 2. Moreover, for 
k > d + 1 there is no d-dimensional lattice that admits 
a A:-simplex state without fine-tuning of the exchange 
couplings. The implication is that for fc > 2 a much 
wider range of ground states are possible in the large- iV 
limit, including spin liquid states. Unless stated other- 
wise, when discussing specific lattices we consider the 
case of nearest-neighbor exchange only. 

We shall first discuss fc-simplex VCS states for the sim- 
pler case Tic = 1, and then generalize to arbitrary Uc- In 
the large- iV limit, by VCS state we mean a saddle point 
where Xrr' is chosen to decompose the lattice into clus- 
ters. Each lattice site belongs to exactly one cluster, 
and any two sites in the same cluster are connected by 
Xrr' 7^ along some path of bonds (they need not be di- 
rectly connected). Each cluster must contain some multi- 
ple of k lattice sites, since otherwise the cluster will not be 
a singlet. In a fc-cluster state, every cluster contains ex- 
actly k sites. A fc-simplex state is a fc-cluster state where, 
within each cluster, each site is (directly) connected to 
every other site by a single bond with J7rr' = Jmax (see 
Fig. [5]). Just as for VBS states, on a given lattice there 




FIG. 2. Illustration of 3-cluster and 3-simplex VCS states on 
the triangular lattice (for ric = 1). Xt-t' is nonzero on the 
highlighted bonds and zero elsewhere. Both states (a) and 
(b) are 3-cluster states. State (b) is a 3-simplex state and is 
a = 00 ground state of the = 3 triangular lattice model. 
State (a) is not a 3-simplex state and therefore has higher 
energy than (b) following the discussion in the text. 

A, A 

\ '/ \ '/ \\ '/ 

A A A 

V " '/ V " '/ 

V '/ V 7 

Xv \/ 

A A 

V '/ 

\v 

FIG. 3. An example fe-simplex state with Uc = 2 and fe = 3, 
on the kagome lattice. Simplices of one color are the triangles 
marked with solid lines (red online), and those of the other 
color are triangles marked with dashed lines (blue online). 
This state was discussed (for = 3) in Ref. [sil . 



can be many different fc-simplex states with the same 
N —i' 00 energy. It is expected that 1/iV corrections will 
select a particular ordered pattern out of this degenerate 
manifold, again precisely as for VBS states ^21, 

To generalize fc-cluster states to nc > 1, we consider 
only diagonal Xr'i-' (^^^ //°^). That is, we consider 

Xt' = S^'^Xrr' (no sum). (83) 
= <5"V" (no sum). (84) 

For each a — 1, . . . , ric, Xrr' chosen to give a fc-cluster 
decomposition of the lattice, resulting in Uc different fc- 
cluster decompositions. A fc-simplex state occurs where 
each fc-cluster decomposition is also a decomposition into 
fc-simplices; an example of a ?ic > 1 fc-simplex states is 
given in Fig. [3] Such states were considered for fc = 2 in 
Ref. and also as exact ground states of special models 
for a variety of ric and N in Ref. 51. 

Focusing on a single color (say, a = 1) and a single 
cluster, and choosing — and xl-r' ^ ~X (for bonds 
within a cluster), the fermionic part of the mean- field 
Hamiltonian in a fc-simplex state is 

^"^'^'^^ = -xE/^i«/wia- (85) 
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The lowest single-particle energy is e£ = — (fc — l)x; the 
k—1 other eigenvalues are degenerate and take the value 
(u = X- The ground state is obtained by filling the lowest 
level in all clusters, and it is easy to see that in this state 
the saddle-point condition {n'f) — mS'^^ is satisfied. This 
state satisfies both the conditions for saturation of the 
lower bound Eq. (|68p . and this is easily verified by direct 
computation of the energy. 

More general /c-cluster states do not saturate the 
bound. To illustrate this, consider for simplicity ti^ = 1 
and a lattice where either Jrri = J^aax, or jTr-r' = 0. 
Consider a /c-cluster state where all the clusters are iden- 
tical and each cluster contains bonds with nonzero 
exchange. It can be shown that the energy of each clus- 
ter is Ec = —NNbJma.^/k'^ (Appendix iPjl . so the total 
energy is then 



MFT 



-Er 



fc3 



(86) 



This attains the bound only if the number of bonds is 
maximum, that is — k{k — l)/2 - but this is precisely 
the condition that each cluster is a fc-simplex. 

As mentioned above, while most lattices admit a VBS 
state, this is not the case for fc-simplex states with k > 2. 
For example, the square and honeycomb lattices admit 
VBS states but no fc-simplex states with fc > 3. The 
triangular (Fig. ^ and kagome lattices admit both VBS 
and 3-simplex states, but lack fc-simplex states for fc > 
4. The three-dimensional pyrochlore lattice of corner- 
sharing tetrahedra admits 4-simplex states, but no fc- 
simplex states for fc > 5. Going beyond specific examples, 
for a d-dimensional lattice, fc-simplex states with fc > d-l- 
1 are impossible, unless the exchange couplings are fine- 
tuned. To see this, consider the fc points of a simplex in d- 
dimensional space. Any pair (r,r') of these points must 
have J^rr' = JTmax! this can only be achieved without fine- 
tuning if space group symmetry forces all the exchange 
couplings to be equal. This can occur only if the points of 
the simplex are mutually equidistant, and there can be at 
most d + 1 mutually equidistant points in d-dimensional 
space. 

While on a given lattice there may be other states 
that saturate the bound even when no fc-simplex VCS 
states exist, for large enough fc saturation is impossible. 
To illustrate this, consider again a lattice where either 
J^rr' = JTmax, Or J7rr' = 0, and let Afb be the total num- 
ber of bonds in the lattice with nonzero exchange. We 
can obtain a lower bound on the energy by treating each 
bond as an isolated system, calculating the resulting two- 
site ground state energy, and summing over bonds. In 
Appendix [C] it is shown that the ground state energy of 
an isolated bond is — nciVj7rnax/fc^, so we have 



UcNNbJn, 
fc2 



(87) 



This bound is more strict than Eq. (|68p when fc > 
2Afb/Ns -f 1, so saturation of Eq. (|68p is impossible for 
such values of fc. 



V. LARGE-iV LIMIT: BIPARTITE LATTICES 
A. Bipartite lower bound 

We now derive a stricter lower bound on the mean- field 
energy that holds for bipartite lattices. As in Sec. IIV A[ 
we consider the mean-field Hamiltonian at general Uc, but 
now on a bipartite lattice. Precisely, we divide the lattice 
into two sublattices A and B of equal size so that t^Trr' 
is only nonzero when r and r' lie in different sublattices. 
We first use the inequality E'mft > -^-mft precisely as in 
Sec. IIV Al The bipartite structure allows us to obtain a 
stricter bound on E^pj,. We recall that 



Ha 



N 



MFT 



\r.r') 



tr(Xr, 



u 



K- 



fr 



(89) 



The crucial observation is that, for a bipartite lattice, Hk 
obeys sublattice symmetry, where TIk — ^ —T~Lk under the 
operation 

fraa V ^ A 
fraa T (E B 

Again we let C be the set of UcmNg occupied levels. 
Now, however, we define the set U to be the image of 
C under the sublattice operation. The set U clearly con- 
tains only empty levels. We denote the set of the remain- 
ing nc{N — 2m)Ns levels by Ai. Levels in Ai are empty 
and have energies intermediate between those in C and 
U. We define averages of and over these sets as 
before. 

As in Sec. IIV Al we have Ej^ — ncmNs[e]£, and we 
need to relate [e]^ to [e^]. We have 

^ 2n 'Tl 2i ™ r 2l (N - 2m) 

2to 2i {N~'2m) 
= £ £. H 



> 



2m , 9, 2m , 



N ' - N 
Since [e]c is negative, this implies 



(90) 
(91) 
(92) 



(93) 



From this point, we can precisely follow the steps of 
Sec. IIV Al to minimize the lower bound on i?MFT- 
this case we obtain the stricter bound 

^ UaNNsJn,^^. (94) 



Emft > 



4fc 



68t when k = 2, and is 



This bound is equivalent to Eq. 
stricter when fc > 2. 



B. Saturation of the bipartite bound 

Here we state the necessary and sufficient conditions 
to saturate the bipartite bound, and give examples of 
fc-cluster VCS states that achieve saturation. 
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(c) 


























(b) 
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FIG. 4. Cluster states {ric = 1) with energies saturating the 
lower bound Eq. (I94[) on the square lattice, for A: = 2 (a), A: = 3 
(b) and k — 4 (c). Xrr' has constant magnitude on the dark 
bonds and is zero on the others. In the fc = 3 state, the flux 
through each six-site plaquette is tt, while it is zero for each 
four-site plaquette in the fc = 4 state. For each value of fc, in 
the N = oo limit, every tiling of the square lattice by the type 
of clusters shown is a ground state. This large degeneracy 
is expected to be lifted upon computing perturbative 1/N 
corrections to the ground state energy.— 



The bound Eq. (IMl) is saturated if and only if each of 
the following two conditions hold: (1) eg is constant over 
each of the sets £ and U, and eg = in 7W. (2) The color 
density n"'' calculated using Hk satisfies the condition 
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tr(/irfir) = 0. 



(95) 



FIG. 5. Illustration of two N — oo cluster ground states on 
the square lattice for = 2 and fc = 4, which saturate the 
lower bound Eq. (|94|) . Square clusters of one color are marked 
with solid lines (red online), while those of the other color are 
marked with dashed lines (blue online). Any configuration 
where clusters of the two colors separately tile the lattice is a 
N — oo ground state - as in the ric = 1 case, the degeneracy 
among these states is expected to be lifted upon computing 
perturbative 1/A'^ corrections to the ground state energy. 



Whenever a given lattice admits a Uc = 1 cluster state 
saturating the bound, it is easy to see that the same 
lattice {i.e. same set of exchange couplings Jrr') also 
admits Uc > 1 cluster states saturating the bound. These 
He > 1 states have diagonal x^r' ^ ^'i- |53)) . and each 
Xrr' chosen to give a cluster decomposition of the type 
that saturates the bound for rif, = 1. Examples of such 
states (for k — 4 and ric — 2) are illustrated for the square 
lattice in Fig. [5] 



VI. LARGE-iV RESULTS ON SQUARE LATTICE 
AND NUMERICAL GROUND STATE SEARCH 



The proof of this statement follows that given for the 
more general bound in Scc. lIVBl As before, condition (2) 
comes from saturation of the inequality -Emft ^ ^mft! 
since nothing in this inequality depends on the bipartite 
structure, the proof of condition (2) is identical to that 
given before. As before, it is trivial to see that E^Y^r^ — 
-Ebound if and only if condition (1) holds. 

As before, saturation of the bipartite bound is impossi- 
ble for large enough k. Again we consider a lattice where 
either Jrr' = Jmax or Jrr' = 0, and let Mb be the total 
number of bonds in the lattice with nonzero exchange. 
For k > 4Afb/Ns, the bound Eq. dHZ]) is stricter than 
Eq. ([M|) . so saturation is impossible for such values of fc. 

Since a flat energy spectrum of the mean-field Hamil- 
tonian is necessary to saturate the bipartite bound, we 
expect that it will only be saturated by VCS states. 
VCS states saturating the bound on the square lattice 
for Tic = 1 are shown in Fig. 2] and were also reported in 
Ref. [3^ For fc = 2 the bound is saturated by any dimer 
state, and for fc = 4 it is saturated by 4-cluster states of 
the type shown. For fc = 3 the bound is actually satu- 
rated by a class of 6-cluster states. 



In this section we focus on the square lattice, and in 
particular on the case fc > 5. The discussion of Sec. IV Bl 
above establishes that, for fc = 2, 3, 4, the large- ground 
states on the square lattice are VCS states, of the type 
shown in Figs. U] and [S] We know of no cluster states that 
can saturate the bound for fc > 5 on the square lattice, 
and we conjecture that saturation is impossible for such 
values of fc. In this situation it is very challenging to 
rigorously determine the large-A^ ground state, a prob- 
lem we do not currently know how to solve. Instead, we 
employ a systematic numerical search for ground states, 
which, while not foolproof, allows us to determine the 
ground state with some confidence. 

Here we first describe our numerical self-consistent 
minimization (SCM) procedure, which we developed and 
employed in Ref. ^ for the case Uc = 1. A very simi- 
lar procedure was later used by M. Foss-Feig and A. M. 
Rey to study the Kondo lattice model, in collaboration 
with one of us (M.H.);^^ and subsequently with both 
of us.i^^ Due to the local constraint, the SCM proce- 
dure is not simply a trivial iteration of a self-consistent 
equation, and to our knowledge it has not been used pre- 
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viously by others; therefore, we shall describe the SCM 
procedure here in some detail. Following this discussion, 
we shall describe the results of SCM on the square lattice 
for ric — 1, 2. 



A. Self-consistent minimization procedure 

We first describe the SCM algorithm in the simpler 
case of He — 1; modifications in the ric = 2 case are 
described below. The basic idea is simply to iterate the 
self-consistency condition Eq. (|13|) . However, if this is all 
one does, then the fermion density will be non-uniform 
and Eq. will be violated. Instead, the idea is to 
iterate Eq. (|13p within a constrained set of Xrr' and fir, 
so that Eq. is always satisfied. To accomplish this, 
the algorithm proceeds as follows: (1) An initial Xrr' is 
chosen randomly. In our calculations, we chose Xrr' = 
\Xrr'\e^'^^^' , where |x| was chosen in the interval [0.03, 
0.18] and </> in the interval [0,27r], both with a uniform 
distribution. (2) Given Xrr', the potential fir is chosen 
so that Eq. (IT4)) is satisfied. We describe below how this 
is done. (3) A new set of x fields is calculated by 



Xr 



N 



-if'', f )■ 



(96) 



(4) We return to step 2, and iterate until the ground state 
energy converges. In practice, we run this procedure for 
500 iterations, by which time the convergence is observed 
to be excellent. 

To improve the efficiency of the algorithm as well as its 
convergence behavior, it is desirable to restrict Xrr' and 
fir to vary within a unit cell, which is then periodically 
repeated to form a larger lattice, with periodic boundary 
conditions. The translation symmetry generated by the 
unit cell primitive vectors allows us to exploit Bloch's 
theorem, further increasing the efficiency. Since different 
unit cells can accommodate different candidate ground 
states, a variety of different cells need to be considered 
separately. 

SCM is indeed a minimization procedure for the 
ground state energy E'mft, as the energy is non- 
increasing for each iteration. To see this, suppose we 
have some Xrr' and fir obtained after step 2. In general 
this is not a saddle point, but Eq. (ITil) is satisfied. We let 
-ffiviFT be the mean-field Hamiltonian defined in terms of 
X and fi. We let xlrr' and fi!r be the fields obtained at 
the next step of the SCM procedure, and -f^MFT i^ 
mean-field Hamiltonian defined in terms of the primed 
fields. We have 



Xrr' 



\Jr\ 



N 



if r' af ra) ' 



(97) 



where the expectation value () is taken in the ground 
state of -ffiviFT- The potential fi'r is chosen so that 
ifrafra)' = whcrc the primed expectation value is 



taken in the ground state of H^pr^. We have 

£^MFT = {HmFt} = N ^'-!-[|Xrr'p-(x;r'Xrr'+C-C 

*Jri~' 



Next, we have the variational upper bound 

-^MFT = (^MFt)' ^ (^MFt) = ~^ 

[r.r') 

Therefore the change in energy satisfies 

^ / \Xrr' — Xrr' P 



\Jrr' 



(98) 



(99) 



(r.r') 



\Jrr' 



< 0; (100) 



that is, the energy is non-increasing for every step of the 
SCM procedure. This means that when the procedure 
converges (in practice it almost always does), it converges 
to a saddle point which is a local minimum of the energy. 
There is no guarantee, however, of a global minimum, so, 
in order to have any confidence that a particular state is 
the global minimum, it is necessary to run the procedure 
many times with different random initial states. 

While the other steps of the algorithm are very simple, 
choosing the potential fir in step 2 requires a more de- 
tailed discussion. The basic idea is to use linear response 
theory to find the change in potential needed to achieve 
a desired change in the fermion density. Going into step 
2, we have fields Xrr' and fir, which can be used to con- 
struct iJMFT- The density will not in general be uniform, 
and we define 



IrO — {ftafra), 



(101) 



where the expectation value is taken using the ground 
state of HyiFT- Suppose the potential is changed by fir — >■ 
fir + Sfir ■ To first order in 6 fir , the change in the density 
is 



Sflr — Xrr' 5 fir 



(102) 



where Xrr' ~ Xr'r = ^rr' i^V definition) the density 
response function evaluated in real space and at zero fre- 
quency, which, using standard results of linear response 
theory, can be calculated from the single-particle wave- 
functions and energies of Hmft ■ While straightforward, 
calculation of Xrr' is the most computationally expen- 
sive step of the algorithm, and must be implemented with 
attention to efficiency. At this point, the idea is to set 
Srir — m — Uro (the deviation between the original and 
desired densities), and invert Eq. (|102p to find Sfir- 

In practice Xrr' is not invertiblc, because the density 
does not change under a uniform shift of fir in the canon- 
ical ensemble. Instead we proceed by diagonalizing Xrr' ■ 



E 



(no sum on a). (103) 
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Here the eigenvalues of Xrr', labeled by a, and 

Ura are the orthonormal eigenvectors. If we expand 
Srir and Sfir in the basis of eigenvectors, we can rewrite 
Eq. (UnH) as 

Sua = XaS^a (no sum On a). (104) 

We invert this by simply ignoring eigenvectors with — 
0, and choosing 

r,, J Sria/Xa , Xa 7^ /'inc;^ 

" \ , = ■ ^^^^> 

This is easily converted back to a result for Siir ■ 

What we have obtained is a linear extrapolation for 
dfir, and the basic idea at this point is to proceed by 
replacing fir fJ-r + ^fJ-r, and iterating the procedure 
until the density is uniform. This is, in fact, just a 
multi-dimensional Newton's method for finding a zero 
of {rir — m) = While such a method has 

good local convergence properties (i.e. starting suffi- 
ciently close to the zero), the global convergence prop- 
erties are poor. However, this can be improved by very 
simple modificationsA22i We define the merit function 
E — X]r("''' ^ rnY , and demand that each change in 
fir decrease £. If fir ^ fJ-r + Sfir actually increases £, 
then we try the smaller step fir ^ fir + \5fir, where 
< A < 1. This is guaranteed to decrease £ for 
sufficiently small A; in practice, we use the sequence 
A = 1,0. 5, 0.4, 0.3, 0.2, 0.1, 0.09, 0.08,..., 0.01, 0.009,..., 
and give up (simply moving on to step 3) after 1000 at- 
tempts. In practice it is only rarely necessary to give up; 
even when it is necessary, step 2 is successful in later it- 
erations, and convergence still occurs. For each iteration 
of step 2, this process of choosing a new 5fir by linear 
extrapolation is repeated 10 times, or until £ < 10^^". 
This tolerance for £ is usually achieved after only a small 
number of iterations, and is virtually always achieved by 
the end of a run (500 iterations). 

Rarely, it happens that £ is of order unity after a sub- 
stantial number of iterations, and the algorithm either 
converges extremely slowly or fails to converge. To avoid 
this problem, when £ > 1 any time after 10 iterations, we 
abort the calculation and start over with a new random 
initial condition. 

We now describe how the SCM procedure is modified 
to handle ric — 2. The initial set of x^r' chosen making 
use of the singular value decomposition 

Here di and d2 are each chosen in the interval [0,0.2] 
with a uniform distribution. U and V are both random 
U(2) matrices, chosen from a uniform distribution on the 
U(2) manifold. 

The algorithm itself proceeds via the same four steps 
outlined above. Only in step 2 are the modifications at 
all nontrivial: We have to choose /i^'' to satisfy {fir^) = 



mS""^. We proceed precisely as above using linear re- 
sponse theory, except that now the necessary linear re- 
sponse equation has a matrix structure and is 

r' c,d 

Since both 5n and Sfi are Hermitian, is is convenient to 
expand them in a basis of Hermitian matrices labeled by 
A, i? = 0,...,3-a convenient basis is the identity matrix 
{A — 0) plus the three Pauli matrices {A — 1, 2, 3). This 
allows one to recast Eq. (|107p in the form 

6n^-Y.T.^rJ^^f^r- (108) 

r' B 

Here, it can be shown that = ^r'r and is real, so 

the response function can be diagonalized as described 
for Tic — 1. Finally, the merit function £ also needs to be 
modified, and we choose 

r r A=l 

B. Results of SCM 

We now describe the results of SCM on the square 
lattice for both Uc = 1,2. The ric = I results were re- 
ported in Reflssl. Following the protocol described below, 
we studied k = 5,6,7,8 for ric = 1, and k = 5,6,7 for 
TT-c = 2. (The numerics become more time-consuming 
with increasing k and ric.) We also checked that SCM 
indeed produces exact ground states (guaranteed by sat- 
uration of lower bounds) for smaller values of k. 

For each value of k and ric noted above, we consid- 
ered all unit cells of rectangular geometry containing fc^ 
or fewer lattice sites, excluding cells of unit width for 
technical reasons. A unit cell of dimensions £x x £y is pe- 
riodically repeated to fill the lattice using Bravais lattice 
vectors R = ixX + iy-y. (Note that other choices of Bra- 
vais lattice vectors are possible ~ we made this restriction 
for the sake of simplicity and limited computation time.) 
The lattice itself has periodic boundary conditions and 
dimensions x Ly. Letting L = mm{Lx, Ly), we always 
considered i > 40 for fc = 5, i > 36 for fc = 6, L > 42 for 
fc = 7, and L > 40 for fc = 8. While in some cases we also 
considered larger system sizes, a more systematic study 
of finite-size effects would be desirable, but we have left 
this for future work. For each unit cell size, we ran the 
SCM procedure 30 times, using a different random initial 
condition each time. 

For ric = 1 and 5 < fc < 8, we found the ACSL to be 
the ground statci^ For ric — 2 and fc = 5, we found the 
ground state to be a rather complicated inhomogeneous 
state that we have not fully characterized. On the other 
hand for fc = 6, 7 we found that the nACSL and dCSL 
are degenerate ground states. 
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VII. DISCUSSION 

We analyzed a variety of SU(-/V) symmetric Heisen- 
berg models in two dimensions on the square lattice and 
gave arguments that topologically ordered spin liquids 
are among their ground states. In view of their potential 
realization with alkaline earth atoms placed on optical 
lattices, we now summarize what we know about realis- 
tically achievable SU(A^) Heisenberg models. Following 
that discussion, we conclude by mentioning some direc- 
tions for future study. 

The Heisenberg models with Uc = 1 can be obtained 
simply as a large-?/ limit (Mott insulator phase) of a Hub- 
bard model representing alkaline earth atoms hopping on 
a lattice with m atoms (in their ground electronic state 
g) per site. Such Heisenberg models are within the reach 
of experiment i^iSi The main issue is temperature, since 
the achieved temperature in experiments is in the range 
f/U < ksT < U, and not ksT < t^/U {t is the Hub- 
bard hopping) necessary for observing effects of magnetic 
exchange. Yet this is similar to the issues encountered 
in studying the SU(2) Hubbard model with cold alkali 
atoms, and currently a significant amount of effort is be- 
ing spent trying to devise techniques to lower the tem- 
perature of Mott insulators. Assuming this is done, the 
study of the Uc = 1 Heisenberg model will be possible in 
the future. 

We summarize what we know about the Uc — 1 Heisen- 
berg model in Fig.|6l On the horizontal axis of this figure, 
we plot m, the number of atoms in the same electronic 
state g per site. On the vertical axis we plot fc, which is 
k = N/m. The dashed-dotted line represents roughly the 
curve km =10. The significance of this curve lies in the 
fact that km — N and iV = 10 is the largest experimen- 
tally achievable N. Therefore, all the points on the plot 
which lie above the curve km = 10 cannot be reached 
experimentally while those below the curve can. The ac- 
tual curve on Fig. [S] is corrected to take into account that 
k and m are integers. 

We emphasize that any iV < 10 is within reach of an 
experiment. Indeed, working with ^'^Sr, for example, one 
can selectively populate its nuclear spin states so that 
only a subset of those are populated with a total number 
of populated states equal to A^.?i At the same time, we 
expect that m — 1 and m ~ 2 columns of the figure are 
easiest to reach, as higher m will likely experience losses 
due to 3-body recombination. 

At m = 1 and k = 2, the ground state is of course the 
Neel state. There is also evidence for magnetic order at 
m = 1, k — 3^- and at m = 1 and k — 4.^- For k — 2 and 
m > 3, it is believed that the ground state is a valence 
bond solid. This is established by quantum Monte Carlo 
for m = 3, 4ji2i and is proven in the limit m — > oo4^ 

In addition to that, in this paper we proved that at 
m — >■ oo, fc < 5, the ground states are valence cluster 
states, of which valence bond solid is a particular exam- 
ple. Finally, we have shown that at fc > 5 and at least for 
fc < 8, and possibly for fc > 8 as well, and at m — (X), the 




1 2 3 4 5 6 m 



FIG. 6. Phase diagram of the SU(A^) Heisenberg model in 
two dimensions on the square lattice with = 1 and with 
— mk. In terms of an underlying Hubbard model, m is 
the number of fermions per site, while k is the inverse filling. 
Regions where there is substantial evidence for a given ground 
state - or where the ground state is known - are shaded. The 
Abelian chiral spin liquid (ACSL) and valence cluster state 
(VCS) regions on the right are established by our large-A^ 
analysis; the boundary between these regions in large-A'^ is 
shown by a dashed line. For = 2, m = 1 the Neel state is the 
well-known ground state. There is also evidence for magnetic 
order at fc = 3, m = 1'^^ and = 4, m = 1.*^ Valence-bond 
solid (VBS) order (which is a type of VCS) was found for 
k = 2 and m = 3, 4:^^ The dashed-dot line separates the 
range of parameters beyond the reach of current experiments 
(above and to the right of the line) and the range within the 
reach of the experiments (below and to the left of the line). 
The experimentally relevant part of the phase diagram with 
the greatest potential for novel ground states - in particular, 
the Abelian chiral spin liquid - is indicated with a question 
mark. 



ground state is the Abelian chiral spin liquid. The rest 
of the phase diagram remains to be filled in. Of course 
other phases not discussed here may well be present, and 
there is some evidence this is the case, in particular at 
k = 2,m = 2M 

The experiments will be conducted at m = 1 or m = 2, 
and at k as large as 10. The ground state of the Heisen- 
berg model under these conditions is not known; this is 
represented by a question mark in Fig. [B] Wc believe it is 
unlikely that the Neel state can survive to large k, even 
at m = 1. Indeed, as discussed earlier, the amount of 
frustration increases with increasing ki^ What happens 
in this region needs to be investigated further. Unfortu- 
nately, numerical study is difficult, especially since these 
models [except when k = 2 (Ref. Il04l )] suffer from the 
quantum Monte Carlo minus sign problem, even on bi- 
partite lattices, in both world-line and fermion determi- 
nantal approaches. However, it may be possible to obtain 
useful information from analytical and density matrix 
renormalization group studies of quasi-one-dimensional 
systems. Ultimately, experiment will need to tell us what 
happens in this part of the phase diagram. An intrigu- 
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ing possibility is that the phase boundary which has be- 
tween k — A and k — 5 extends all the way from large 
m to ni — 1, thus making the experimentally accessible 
m — \, k > A regime a chiral spin liquid. 

We note that, while we only considered integer k, some 
non-integer values of k are possible. For example, m — 2 
and N = 5 corresponds to k — 5/2, and a well-defined 
large- limit with k = 5/2 certainly exists. We did not 
consider such values of k first for simplicity, and second 
because non-integer k requires m > 2, making experi- 
mental accessibility somewhat less favorable. Nonethe- 
less, it would be interesting to study the large-A^ limit 
for non-integer values of k in future work. 

A similar phase diagram can be discussed at ric = 2 
where the Abelian chiral spin liquid will be replaced by 
the non- Abelian chiral spin liquid (or by the doubled chi- 
ral spin liquid). 

Supposing that some of the topological liquids dis- 
cussed here do indeed occur for physically realizable 
SU(iV) spin models, it will be an interesting question how 
to actually observe fractional or non- Abelian statistics in 
these systems. This is especially so given the intense 
interest in topological quantum computation using non- 
Abelian particles. We expect that holes, the excitations 
obtained by removing an atom from the system, should 
split into spinous and holons. The holons may be local- 
ized near a given site by an external potential and, at the 
same time, they obey fractional or non- Abelian statistics 
depending on which topological liquid we are considering. 
Therefore braiding may be achieved by manipulating the 
holons via the external potential, and this is a route by 
which fractional and non- Abelian statistics may be ob- 
served. While some further details along these lines are 
given in Appendix [Fl many questions remain open, and 
we feel this constitutes an interesting direction for future 
work. 

Other directions for future study include investigation 
of the projected wavefunctions for the various topolog- 
ical liquids, which we discussed only briefiy in Sec. IIIII 
Given the difficulty of unbiased numerical study in these 
systems, such wavefunctions may be useful as variational 
states to gain understanding of the phase diagram away 
from the large- iV limit. Finally, another potentially inter- 
esting problem is a careful study of the dCSL edge states, 
which may be topologically protected as mentioned in 

Sec. Unci 
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Appendix A: Alkaline earth atom Hubbard and spin 
models 



Here we briefiy review the Hubbard model describing 
fcrmionic alkaline earth atoms in optical lattices. We 
focus on two kinds of Mott insulating states, in which 
the spin models we study are the simplest description 
capturing the essential physics. A more extensive and 
detailed discussion of fermionic AEA in optical lattices, 
and the rich variety of strong correlation physics that can 
be realized in these systems, can be found in Ref. [2ll 

A single alkaline earth atom has a ^So electronic 
ground state. (Recall that the subscript on the right 
is J, the electronic angular momentum, so this state has 
J ~ 0.) The nuclear spin can be as large as / = 9/2 in the 
case of ^^Sr. Other important examples are -'^^^Yb and 
"^Yb, with / = 1/2 and / = 5/2, respectively. While Yb 
is not an alkaline earth, it has the same configuration of 
outer electrons, and all the discussion here applies equally 
to alkaline earths and to Yb. Also important for our 
purposes is the "^Pq lowest electronic excited state, which 
has a very long lifetime on the order of 100 s. These two 
electronic states can be subjected to optical lattices of 
different strengthJ^ 

Interactions between two atoms in any combination of 
these electronic states, which arise from collisions in the 
s-wave channel, are expected to respect a large SU(A^) 
spin rotation symmetry, where A'^ = 2/ -I- 1 is the number 
of nuclear spin levels per atomi2ii2^ The symmetry arises 
because such atoms have J = 0, and due to the resulting 
quenching of hyperfine coupling, the nuclear spin is essen- 
tially a spectator in the collision between two atoms, and 
only participates via Fermi statistics. The SU(iV) sym- 
metry is not exact but is expected to hold to an excellent 
approximation. A rough estimate is that, for two ground 
state atoms, SU(A^)-breaking effects are 10^^ times the 
strength of the SU(7V)-symmetric interaction.— For two 
excited state atoms, the strength of SU(7V)-breaking is 
estimated to be 10"'^. 

We now suppose that the atoms are subjected to an op- 
tical lattice potential deep enough that a description in 
terms of a one-band Hubbard model is appropriate. We 
introduce creation operators c^g^ and cl^^ for the ground 
state {g) and excited state (e) atoms, respectively. Here, 
r labels the lattice site, and a = 1, . . . ,N labels the z- 
component of nuclear spin. (We shall find this notation 
more convenient than = —I, ...,/, due to the SU(Af) 
symmetry.) To describe the system we consider the most 
general Hubbard model with SU(Af) symmetry, nearest- 
neighbor hopping, and on-site interactions. It is also im- 
portant to note that the numbers of ground state and 
excited state fermions are separately conserved, due to 
the long lifetime (treated here as infinite) of the excited 
state fermions, and energy conservation. The Hamilto- 
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man is 

H = -tgY^ {clg^C^,g^ + H.C.) -t^Y^ (4ea Vea + H-C.) 
{rr') {rr') 

+ 51 (^"-rg + ^^le + UegUrgrire) 
r 

-■^e,E^^^(r)53Jr). (Al) 

r 

The sums in the first two terms are over nearest-neighbor 
bonds. We have introduced the fohowing number and 
spin operators: 

a 

with corresponding expressions for rire and S^p{r). The 
on-site interaction parameters C/^g, Uee-, Ueg and Jeg are 
proportional to hnear combinations of the four indepen- 
dent s-wave scattering lengths characterizing collisions 
among the atoms.— 

The SU(A'^) spin symmetry acts on the fermions as 
follows: 

4ea ^ C^a/JCL/S- (A4) 

Here, U is an arbitrary SU(A^) matrix. The fermions thus 
transform in the fundamental representation of SU(A^). 

We shall consider Ugg > 0, which is known to be the 
case for *^Sr and ^^'^Yb. In both cases the corresponding 
scattering length is about 100 ao, which corresponds to 
rather large repulsive interactions . The sign of the 
interspecies exchange interaction Jeg is not yet known 
and may be either ferromagnetic (positive) or antifer- 
romagnetic (negative); this is likely to depend on the 
atomic species. If one ground state atom and one excited 
state atom share the same site, antiferromagnetic (ferro- 
magnetic) Jeg favors antisymmetric (symmetric) combi- 
nations of their nuclear spins. 

We consider two types of Mott insulators. The simpler 
of the two is realized using only ground state atoms, at 
an integer filling of m atoms per site. While m — 1 best 
avoids issues of three-body loss, we consider general m 
because it is needed for the large- limit. In this case, 
the Hubbard model contains only the tg and Ugg terms, 
and when tg <^ Ugg the standard degenerate perturba- 
tion theoryiS? gives the spin model 

Hspin ^JJ2 Sc.0ir)Si3c.{r'), (A5) 

{rr') 

where J — 2t'^/Ugg, and we have defined 

S^pir)^Si^{r) + S^^p{r). (A6) 

In this case, since no excited state atoms are present, 
Sapi^) = S^^Jr). The spin at each site transforms in 



the m X 1 irreducible representation of SU(Af) (Fig. [Ij; 
this simply expresses the fact that the nuclear spins of 
the identical fermions are combined antisymmetrically. 
For m = 1, the spin transforms in the fundamental rep- 
resentation of SU(A^), and when N = 2 this is simply a 
5 = 1/2 spin. 

The second type of Mott insulator is related to 5* = 1 
Mott insulators of SU(2) spins. It is realized with m 
ground state atoms and m excited state atoms on each 
site. We consider Jeg > 0, in which case the single-site 
ground state is associated with a m x 2 tableau. This 
can be seen by first viewing the ground state atoms as 
forming a spin transforming in the to x 1 representation, 
which is required simply by Fermi statistics, and simi- 
larly for the excited state atoms. These two spins are 
then coupled by the Jeg exchange term, and formally we 
need to solve a two-site problem, which is done in Ap- 
pendix [C] We do not consider antiferromagnetic inter- 
species exchange, because this gives a single-site ground 
state with a 2to x 1 tableau, which is of the same type 
obtained with only ground state atoms 

In the simple case to = 1, the single-site ground states 
are of the form 

\^r) = ictgJ,e0+4-gpcleam; (A7) 

that is, the nuclear spins of the two fermions are com- 
bined symmetrically. When N = 2 and to = 1, this 
is simply a 5 = 1 spin. More generally the single-site 
ground states can be obtained from the highest-weight 
state 

\tpj. ) — clgiclei ■ ■ ■ clg,j^cle.fn\0} J (A8) 

where all other single-site ground states can be obtained 
by repeated action on IV'J.™) with appropriate compo- 
nents of Sapir). That is, they are linear combinations of 
states of the form Sap{r)\^^;.'^) , Sapir)S^s{r)\^r^) , and 
so on. 

Again, to = 1 best avoids issues of three-body loss, 
but we shall consider general to. Another potentially 
important loss mechanism is inelastic losses in collisions 
between two excited state atoms. This can be minimized 
by making the lattice for the excited state atoms very 
deep, effectively setting te = 0- 

In Sec. ini the type of SU(A^) spin is specified by the 
two local constraint equations Eq. ([2]) and Eq. ([3]) , and in 
Appendix |B] it is shown that these two constraints imply 
that the spin transforms in the m x Uc representation. 
To make contact with that discussion, we now show that 
single-site ground states of the present Hubbard model, 
transforming in the to x 2 representation, satisfy the con- 
straint Eq. that is 

t;iv.) = o, (A9) 

where 

K - lci,^<bC,ba- (AlO) 
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Here, a,b = e, g, and we formally consider the e, g labels 
as an index transforming in an SU(2) "orbital" space. 
Moreover, cr* are the 2x2 Pauli matrices {i — 1,2,3), 
and \ipr) is a single-site ground state for the site r. [Since 
there are 2m fermions on each site, the constraint Eq. ([2]) 
is obviously satisfied.] 

The constraint Eq. (jA9[) is obviously satisfied for the 
m = 1 state given in Eq. (IA7|) : the wavefunction is an- 
tisymmetric under interchange e -H- g and is thus an or- 
bital singlet. The same holds for the highest-weight state 
iV'r^), since it is built as a product of orbital singlets 
4ga4ea ("^o sum On a). Becausc [S'a^(r),r^] = 0, this 



immediately implies that Eq. (IA9p holds for all single-site 
ground states. 

When te = and tg ^ Ugg , Jeg , the spin Hamil- 
tonian is given by the same form Eq. (jASp . only now 
J = ig/[2([/gg -I- Jeg)]- The degenerate perturbation the- 
ory calculation needed to establish this, unlike in the case 
of only ground state atoms, is not simply a trivial general- 
ization of the familiar calculation for the S — 1/2, SU(2) 
Hubbard model. While the end result of this calculation 
appeared in Ref. [U the details were not presented, so 
we now present them here. 

First we consider a single lattice site, and note that 
the energy of 1^/"^^) (neglecting hopping) is 
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1 



^0 = -jUgg-m^ + -Uee-m + UegTn - Jeg-m. (All) 

By SU(A^) symmetry, this holds for any single-site ground 
state. Moreover, we note that 
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C W^"") 
^rea^rgaWr I 



0, 



(A12) 



which also holds for any single-site ground state by 
SU(A'^) symmetry. 

Now we consider second-order degenerate perturbation 
theory for a two-site problem with adjacent lattice sites 
ri and r2. (In second-order perturbation theory, there is 
no need to consider more than two sites.) We construct 
the effective Hamiltonian by building up its action on 
an arbitrary state |'0ri)l'0r2) low-energy manifold 

(that is, IV'ri) 9.nd IV'ra) arbitrary single-site ground 
states). The energy of the initial state is 2Eq. The in- 
termediate state is obtained by hopping a single ground 
state fermion from ri to r2 or vice versa. We consider 
hopping from ri to r2 so the intermediate state is 



IV'int) 



I'?^ria)l0r2a)' 



where 



(A13) 

(A14) 
(A15) 



Acting on the intermediate state with the on-site part 
of the Hamiltonian and using the identity Eq. (jAl2p to 
evaluate the action of the J^g exchange term, the energy 
of the intermediate state is found to be 



E'int = 2£'o + [/„„ + Je, 



(A16) 



Since the intermediate state is an eigenstate with en- 
ergy independent of the initial state, the effective Hamil- 
tonian is 
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Ua 



Je 



'gg r ^eg 

+ (ri o r2) 



'^'-rigQ'^r2ga(l '^^'^r^gfi'-'rigfi 



V 



(A17) 



where V is the usual projector onto the ground state 
manifold, and the second term in the square brackets 
accounts for the process where a fermion first hops from 
T2 to ri . Because a single hopping process always leaves 
the ground state manifold (since it changes the fermion 
number on each site), we can drop the (1 — P) factor and 
write 



-ti 



gg 



Jeg 



'Pcl^gaC^^gc.cl^gd^r.gd'P + i^l ^ 



^^g 



Je 



(A18) 



where we dropped an additive constant in going to the 
second line. Now, 

SUir) = \Se.p{r) + \[Slp{r) - ^^^(r)], (A19) 

where the second term transforms as a triplet in the or- 
bital space. But the projector V forces every lattice site 
to be an orbital singlet, and therefore 

'PSi0{r^)Sl^{r2)V = \Se,0{r^)Sp^{r2), (A20) 



so that 



H, 



cff 



''g 



2{Ugg + Jeg) 

the result we claimed above. 



Sa.fs{ri)Si:sa{r2), 



(A21) 



Appendix B: Determining the irreducible 
representation of SU(A'^) spin from local constraints 

Here, we show that the local constraints Eqs. (|2l3p im- 
ply that each spin transforms in the m x He irreducible 
representation of SU(Af). Another way to state this fact 
is that the Hilbert space of a single lattice site, sub- 
ject to the local constraints, transforms irreducibly under 
SU(A^) in the m x Uc representation. 

To see this, it is helpful to think of spin and color rota- 
tions as a subgroup SU(iV) x SU(nc) C SU(nciV), where 
the fermions transform in the fundamental of S\J (ucN). 
By fermion antisymmetry, the first constraint [Eq. ^] 
implies that each site transforms in the ricm x 1 repre- 
sentation of SU(nciV). 

To understand the role of the second constraint 
[Eq. Q], we need to understand how this representation 
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FIG. 7. Illustration of Young tableaux occurring in the de- 
composition of the ricm x 1 representation of S\J{ncN) into 
irreducible representations of the SU(A'^) x SU(nc) subgroup, 
for the case ric = m — 2. If, as described in the text, 
we project out the SU(nc) irreducible representation corre- 
sponding to the tableau on the left, then the corresponding 
SU(A'^) tableau is as shown on the right. Note that the rows 
(columns) of the SU(nc) tableau become the columns (rows) 
of the SU(iV) tableau. 



is given by interchanging the role of rows and columns 
of the SU(nc) tableau - see Fig. [7] for an example that 
clarifies the meaning of this statement. 

The constraint Eq. ^ dictates that we keep only the 
terms in the decomposition where ?'gu(„ ) is the singlet 
representation Osu(nc)- Since we have to form the cor- 
responding tableau using ricm boxes, the only possible 
SU(nc) tableau is ric x m, and the above discussion im- 
plies that the corresponding SU(A^) tableau is m x ric- 
It can be seen by directly constructing a highest weight 
state that the representation (mx nc)su(Af)'XiOsu(ne) only 
occurs once in the decomposition. Therefore the con- 
straint gives 

{ficfn X l)su(ncAr) -5> (to X nc)su(w) Osu(„,), (B4) 
the desired result. 



decomposes into irreducible representations of SU(A^) x 
SU(nc)- The decomposition has the general form 



(ncm X l)su(neW) = ^^su{n) 



(Bl) 



This equation expresses the fact that the {ricm x 1) repre- 
sentation of SU(nciV) is a direct sum of irreducible repre- 
sentations of SU(A^) X SU(nc), labeled by i. We will first 
show that, for each term in this decomposition. 



uniquely determines rl 



' SV{N) 



SU(ne) 



and vice versa. 



Focusing on a single lattice site and dropping the site 
label for fermion operators, we consider the following 
(overcomplete) basis states for the {ricm x l)su(racAf) ^^P" 
resentation: 



|ai, ai; 



\ ft ft 



JO). (B2) 



If P{i) is a permutation of the integers i — 1, 
then fermion antisymmetry implies 

|ap(i), ap(i); ■ • ■ ; ap(ricm)7 ap(riem)) 



sgnP|ai,ai;...;a„ 



1 1 



, , ricm, 



(B3) 



where sgnP is the sign of the permutation. Suppose we 
want to project out a particular representation of SU(nc). 
We do this by forming a SU(nc) Young tableau with ricm 
boxes, and associating each box with a color index a^. We 
then follow the usual procedure of first antisymmetrizing 
the ai indices occupying the same column, and second 
symmetrizing those occupying the same row. Because of 
the overall antisymmetry expressed in Eq. (jB3p , when in 
the first step we antisymmetrize the ai indices in a given 
column, we also simultaneously symmetrize the corre- 
sponding set of ai indices. Similarly, the second step an- 
tisymmetrizes those ai indices corresponding to a given 
row. This means that, in the process of projecting out 
a given desired SU(nc) representation, we have also au- 
tomatically projected out a corresponding given SU(A^) 
representation. The tableau of the SU(7V) representation 



Appendix C: Exact ground state energy of two-site 
problem 

Here we consider a problem of two spins at ri and r2 , 
coupled by the Hamiltonian Eq. ([S]). We write J^ra = 
J/N, so that 



n 



J 

N 



Sal3{ri)Spa{r2)- 



(CI) 



We shall calculate the exact {i.e. not large- A'^) ground 
state energy for arbitrary N, m — N/k and ric. 

It is convenient to define the Hermitian spin operators 



r Jraa aBJ'i 



(C2) 



where A=l,..., N'^-l labels the S\J{N) generators T-^. 
These are chosen to satisfy the orthonormality condition 



AB 



tr(r'^T") = -5 
and can be shown to satisfy the identity 
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(C3) 
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Equation (|C4[) can be used to show 
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Now, {f^f = fAfA jg ^i^g quadratic Casimir of 
SU(iV). In a given irreducible representation r this op- 
erator is proportional to the identity, and its eigenvalue 
6*2 (r) can be computed from the structure of the Y oung 
tableau using a formula given in Ch. 19 of Ref. Il09l 
which we now reproduce. Suppose the Young tableau 
has rirow rows, each with length hi {i = 1, . . . , rirow) and 
ricoi columns, each with length ai {i = 1,. . . ,ncoi), and 
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a total of i boxes. Then the eigenvalue of the Casimir is 
given by 



C2{r) 



(C6) 



1=1 



i=l 



Since each spin transforms in the mxric representation, 
we can use Eq. ^ to evaluate (f^)^ = (^4)^- More- 
over, by examining the Young tableaux appearing in the 
tensor product {m x ric) (g) (m x ric), and using Eq. (|C6|) 
to evaluate (T^ +T^)'^ for each tableau, we find that the 
two-spin ground state is the 2m x ric tableau, and that 
the corresponding ground state energy is 



ricNJ 
fc2 • 



(C7) 



We know this must be the ground state energy of the iso- 
lated fc-cluster because it saturates the bound Eq. (|87|) 
provided by the ground state energy of the two-site prob- 
lem. 

We note that this result also holds at any finite N . 
Schematically, this can be seen by noting that the ground 
state is the unique singlet that can be formed from the 
k spins, which can be thought of as a x nc Young 
tableau, which is obtained by vertically stacking the m x 
Uc tableaux for each of the k sites. Any pair of spins can 
then be seen to transform in the 2m x Uc representation, 
which, by the discussion of Appendix [Cl implies that the 
two-site Hamiltonian on the link connecting those sites is 
in its ground state. So the ground state energy is just the 
sum of the two-site ground state energies for each link in 
the graph, which again gives Eq. (|D6p . 



Appendix D: Energy of fc-cluster states 

Here we compute the large- ground state energy of 
a single isolated fc-cluster, a result which is used in the 
discussion of Sec. lIVCl We consider a spin model defined 
on an arbitrary connected graph with k sites labeled by 
s, and with links labeled by l. The exchange energy is 
taken to be equal on all links and is J = J /N . The 
mean-field Hamiltonian is 

i^MFT = ^ 51 tr(xlx^) + tr(A's) + Hf, (D1) 

^ I s 

where T-Lf — T~Lk -\-T~Lv i and the latter two operators are 
constructed as in Eqs. (|19l20p . 
We consider the following ansatz: 

Xf = -S'^'x (D2) 
fif = -S-'zsX- (D3) 

Here, Zs is the coordination number of the site s. We shall 
see that x > upon minimizing the energy with respect 
to X- With this choice, fixing the color and spin quantum 
numbers, the one-particle Hamiltonian that can be read 
off from 7^j7 is proportional to the Laplacian matrix of 
the graph (with positive coefficient). Therefore the single 
particle ground state (for fixed color and spin) has zero 
energy, is unique, and its wavefunction is a constant. The 
unique many-particle ground state oi'Hp is obtained by 
filling this state with kmric — ricN fermions, one in each 
of the UcN possible combinations of color and spin states. 
The mean-field energy is therefore given entirely by the 
constant terms in iJMFT, and is 



Emft = :f — X ~ mricX z 

= — J — X ~ 2mncNbX, 



(D4) 



(D5) 



where Nh is the number of links in the graph. Minimizing 
with respect to x, we find 



-Emft 



fc2 • 



Appendix E: Chiral spin liquid in large-fc limit 

Our demonstration that constant magnetic field with 
a fiux of 27r/fc per plaquette is the lowest energy solution 
to the saddle point equations - on the square lattice, for 
Uc — 1 and 5 < fc < 8 - is purely numerical. A natural 
question which arises in this context is whether this can 
be supplemented by additional analytical analysis. 

The solution to the problem of a particle hopping 
on a square lattice in a constant magnetic field of fiux 
27r/fc cannot be found analytically. Yet it is well known 
that the spectrum consists of fc bands (Landau levels).-^ 
Since the fermions we work with are at a filling fraction 
1/fc (filling all the bands would correspond to N parti- 
cles per site, while we have instead m = N/k particles 
per site), they fill precisely one lowest Landau level. Yet 
the energy of a filled Landau level is not known analyti- 
cally, except at a very large fc where the problem becomes 
effectively continuous. 

Therefore, let us calculate the energetics of a saddle 
point solution with a flux of 27r/fc per plaquette (corre- 
sponding to ACSL) in the limit of very large k and com- 
pare it with other possible states at this fc. Our analysis 
will be for the case Uc — 1, but also applies immediately 
to 71c = 2, since any ric — 1 saddle point can be extended 
to a ric — 2 saddle point of the diagonal form 



A.rr 



S'^^Xrr' (no sum) 



= 5°Vr (no sum). 



(El) 
(E2) 



(D6) 



where each pair {Xrr''f^r) is a nc = 1 saddle point solu- 
tion. The energy is simply a sum of energies of the Uc — 1 
saddle points. We can obtain the nACSL and dCSL sad- 
dle points in this fashion from the ACSL saddle point, by 
choosing x^ and to have the same or opposite mag- 
netic fields, respectively. The other ric — 1 states we 
consider can similarly be straightforwardly extended to 
Uc = 2 states in this fashion. 

We start by choosing the hoppings Xrr' according to 
Eq. First of all, the first term in the mean-field 
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energy is easy to calculate 



{rr' ) 



J 



(E3) 



Here, as before, Ng is the total number of sites on the 
lattice, and 2Ns is the total number of bonds. 

Now let us find the energy of a fully filled Landau level. 
A particle hopping on a lattice with a hopping strength 
X without the magnetic field has the spectrum 



-2xcos(fc^) - 
-4x + X {kl 



2x cos{ky 



k 



'■)- 2i. (hi 



(E4) 



where a small k^, ky expansion was performed (lattice 
spacing is taken to be unity). Looking at the quadratic 
term, we read off the effective mass of the particle m* — 
f/(2x). This gives the cyclotron frequency 



UJ 



B 

m* 



Awx 
k ' 



(E5) 



since the magnetic field is B 
lowest Landau level is then 



27r/fc. The energy of the 



El 



-Ax+^uj 



'4x 



27rx 



(E6) 



For what follows we would like to also calculate the 
1/k^ correction to this result. The corrections come from 
the quartic term in the dispersion, which takes into ac- 
count the deviation of the lattice from the continuum 
limit. The correction to the Hamiltonian describing the 
motion of a particle in a magnetic field in the continuum 
due to this term in the dispersion can be found by min- 
imal subtraction (for example, in Landau gauge), and 
gives 



V 



x_ 

12 



I — 

dx 



2^y\ 
k ) 



dyi 



(E7) 



Considering this a perturbation, the unperturbed wave 
function is given byiii 



V'(fc.) - ( ^ 



^) ). (E8) 



Calculating the matrix element {il^{kx) \ V IV-'(fcx)) we find 



El = -4x + 



2nx 



2fc2 ■ 



(E9) 



This is the energy of the lowest Landau level in the ap- 
proximation up to terms l/fc^. Notice that the Landau 
level remains flat, that is, k^ independent. It is easy to 
see that it will remain flat up to arbitrary order in 1/k. 
This means that the broadening of the Landau level is 
exponentially small in 1/k and can be ignored for the 
purposes of this calculation. 



The total number of particles fllling the Landau level 
is NNs/k, so the mean field energy becomes 



E, 



MFT 



2NsNx^ NNsEl 
J k 

Minimizing this with respect to x we find 



E, 



MFT 



JNN, 



2-K 

T 



TT 



(ElO) 



(Ell) 



Now let us consider alternative states. One alternative 
state is a Fermi surface state, where all hoppings are real 
and equal to x- The energy of such a state is straight- 
forward to calculate. We take particles moving with the 
dispersion given by Eq. (jE4l) . fill all the states at an ap- 
propriate density up to Fermi energy, to find the total 
energy per particle to be 



Ef 



'4x 



k 



3fc2 ■ 



(E12) 



This energy is slightly higher than the energy of the Lan- 
dau level given in Eq. (jE6p . Therefore the energy after 
minimization with respect to x is also slightly higher 



E 



JNNs 



MFT 



2 27r 57r^ 
k 6fc2 



(E13) 



Clearly, Eq. (|Elip is greater than Eq. (|E13p , so the state 
with the uniform magnetic field wins. 

A second alternative state one might consider is a VCS 
state. Suppose the lattice is covered by clusters of exactly 
fc-sites each, each containing Nf, bonds. Within each clus- 
ter Xrr' are constant and equal to x, which is real, and 
Xrr' for bonds connecting the clusters are zero. In Ap- 
pendix |D] it is found that the energy of a single cluster is 
given by Eq. ID6I) . Since the cluster energies simply add, 
and the number of clusters is Ng/k, the total energy is 



i?MFT 

Now we can define N^e by 



JNNgNb 

fc3 



= 2Ng - Nbe, 



(E14) 



(E15) 



so that Nbe is the total number of bonds not contained 
inside some cluster. We have 



Emft — —J- 



.2NNs 



fc2 



J 



NNbe 
fc2 ■ 



(E16) 



Nbe scales with the total perimeter of all clusters. Since 
the perimeter of a single large cluster goes like Vk, Nbe 
is proportional to \/k times the total number of clusters 
Ns/k, or 



'Vk' 



(E17) 



24 



where c is some constant. This gives 



Emft 



(E18) 



Comparing this with the uniform magnetic mean-field 
energy Eq. (lElip as weh as the uniform hopping given 
by Eq. (jE13p . we see that the magnetic field mean- field 
energy is again the lowest at large k. 

The arguments presented here do not prove that the 
uniform magnetic field is the lowest energy solution. 
That is demonstrated instead by the numerical solution 
of the mean-field equations. However, they do give a feel 
and perhaps some intuition as to why this solution wins 
over some of the possible alternatives. They also sup- 
port the idea that chiral spin liquids are good ground 
states not just a few intermediate values of fc, but also 
for larger k. Therefore we conjecture that the ACSL is 
the large- TV ground state for ric = 1 and all fc > 5, and 
that the nACSL and dCSL are the degenerate large- 
ground states for Uc = 2 and all fc > 6. 



Appendix F: Localization and braiding of fractional 
and non-Abelian particles 

One of the most striking properties of the topological 
liquid states discussed in this paper is the presence of 
particles with fractional and non-Abelian statistics. It 
is therefore interesting to discuss how, in principle, such 
particles may be localized and braided, especially in view 
of the intense interest in topological quantum computa- 
tion using non-Abelian particles. Our intent here is not 
to develop a detailed and realistically achievable proposal 
to carry out such a braiding experiment in a cold atom 
system, but simply to discuss in principle how such braid- 
ing may be achieved, and point out some of the issues 
that arise. Development of more detailed proposals is 
an interesting subject for future work. It would be also 
interesting if our discussion can be sharpened by appro- 
priate calculations. For ease of presentation, we focus on 
the case ric — m — 1; generalization to other cases is 
straightforward. 

We shall discuss how one may localize a particle called 
a holon that is spinless but carries the conserved atom 
number. The reason we consider holons rather than 
spinous, is that holons may be localized simply by modi- 
fying the optically generated single particle potential for 
the atoms. To do this, we need to go beyond the Heisen- 
berg spin model, and for greatest simplicity we consider a 
t-J model where strong correlation restricts the number 
of atoms per site to be less than or equal to one. The 
Hamiltonian is 

HtJ ^-tgY.'^ + H.C.) V 



{rr' ) 

J t t 

^ / ^ ^r' a^'^' l^^r(i^'^^'> 
{rr') 



(Fl) 



where c|.„ creates a ground state atom in spin state a on 
site r, and V \s & projector onto the subspace with one 
or fewer atoms on each site. Jrr' bas been replaced by 
J on every bond, and the sum in the first term is over 
nearest-neighbor bonds of the square lattice. Cra is said 
to insert a hole (with spin a) at site r. When there are 
no holes present, this model reduces to the Heisenberg 
spin model with ric = to = 1. Below, we rely on the 
approach of Lee and Nagaosa to discuss this modelJ^ 

To make contact with the description of the topological 
liquid states, we decompose the hole insertion operator 
as 



ra. " Jra^ri 



(F2) 



where f^^ creates a spinon and h\. is a bosonic creation 
operator creating a holon. Spinon and holon densities 
obey the local constraint 



P f + feU = 1 

J ra. J ra ' r 



(F3) 



Assuming the system (without holes) has an ACSL 
ground state, the spinous are low-energy quasiparticles, 
which couple to the Chern-Simons gauge field and thus 
acquire fractional statistics. The holon also carries gauge 
charge, and thus also acquires fractional statistics. It 
should be noted that, in the equations above the holon 
and spinon are formal objects used to microscopically 
represent the t-J model, and these formal objects are not 
the same as the low-energy quasiparticle degrees of free- 
dom. Holons and spinous emerge as low-energy degrees 
of freedom when we study the t-J model starting from 
an appropriate mean-field theory^ii^ and then including 
fluctuations. Since the discussion here is only qualitative, 
and since the needed mean-field theory is very closely re- 
lated to that introduced in Sec.|TTl we shall not introduce 
it here. It suffices to note that, at the mean-field level, 
both spinous and holons are free particles, which are min- 
imally coupled to the fluctuating gauge field upon going 
beyond mean-field theory. 

We now consider introducing the external potential 



SHu = -J2u{r)cl 



(F4) 



and adding a single hole into the system. The sign in 
Eq. (|F4p is chosen so that a negative U{r) is an attractive 
potential for the added hole. Up to an additive constant, 
we may use Eq. (|F3p to re-express the potential as 



5Hu = J2uir)bX. 



(F5) 



We could have equally chosen the potential to couple 
to the spinons and not the holons; the above choice is 
convenient, but is purely a convention. For example, at 
the mean-field level, a change in the saddle point value 
of the Lagrange multiplier field enforcing the local con- 
straint, will apportion the effect of U{r) between holons 



25 



and spinons. Therefore the system dynamically deter- 
mines the effect of the physical external potential U{r) 
on holons and spinons. 

When the hole is added, Eq. IF2I tells us that we both 
add a holon and remove one spinon. (The removed spinon 
should really be called a spinon hole, but for ease of dis- 
cussion we will simply call it a spinon.) We suppose that 
U{r) is negative, appreciable only in a small spatial re- 
gion, and just strong enough to bind a particle. Because 
U{r) couples to the conserved density, we expect it to 
bind a particle carrying atom number —1, but it is not 
obvious whether this particle will be a hole or a holon. To 
understand this, the added holon and spinon will interact 
via some short-ranged potential. This potential may be 
attractive or repulsive. If the holon-spinon potential is 
attractive enough, the holon will be bound to the spinon, 
and they will be localized together by the external po- 
tential U{r). In this case we have localized a hole, which 
is not a fractional particle. On the other hand, if the 
holon-spinon potential is repulsive enough, a holon will 
be localized. In the latter case, we can then manipulate 
the fractional holon by adiabatically changing the exter- 
nal potential U{r). Multiple holons could be created by 



choosing U{r) to be a sum of several localized potentials. 

Since the goal is to create and manipulate a fractional 
particle, what should be done if a hole is localized by 
the external potential? One solution is to apply a time- 
varying Zeeman magnetic field, which will couple to the 
localized spinon and can be used to excite it to a delocal- 
ized state, leaving behind a localized holon. If we do this 
to create a state with several localized holons, the delo- 
calized spinon excitations will induce some errors when 
the holons are braided. However, these spinon excita- 
tions can be made to relax by whatever cooling mecha- 
nism was used to prepare the state of several localized 
holes in the first place. (Finding a cooling mechanism 
capable of achieving this for cold atom Mott insulators is 
a significant unsolved problem. Solving it is a prerequi- 
site for any experiment probing fractional or non-Abelian 
statistics in such systems, which would have to be car- 
ried out at temperatures well below the bulk gap.) While 
some spinons may relax back into the localized states and 
re-form holon-spinon bound states, because these states 
are localized, the rates for other relaxation processes (for 
instance, relaxation into low-energy edge excitations) are 
expected to dominate. 
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